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ABSTRACT 


A  new  approach  is  proposed  for  analysing  the  compressible 
turbulent  boundary  with  arbitrary  pressure  gradient.  The  new 
theory  generalises  an  incompressible  study  by  the  first  author 
to  account  for  variations  in  wall  tesperature  and  freest ream  Mach 
number  and  temperature.  By  properly  handling  the  law -of -the -wall 
in  the  Integration  of  momentum  and  continuity  across  the  boundary 
layer,  one  may  obtain  a  single  ordinary  differential  equation  for 
skin  friction  devoid  of  integral  thicknesses  and  shape  factors. 

The  new  differential  equation  is  analysed  for  various  cases. 
For  flat  plate  flow,  a  new  relation  is  derived  which  is  the  most 
accurate  of  all  known  theories  for  adiabatic  flow  and  reasonably 
good  (fourth  place)  for  flow  with  heat  transfer.  For  flow  with 
strong  adverse  and  favorable  pressure  gradients,  the  new  theory 
is  in  excellent  agreement  with  experiment,  possibly  the  most 
accurate  of  any  known  theory,  although  the  data  are  too  sparse 
to  drew  this  conclusion.  The  new  theory  also  contains  an  explicit 
criterion  for  boundary  layer  flow  separation.  Also,  it  appears 
to  be  the  siaplest  by  far  of  any  compressible  boundary  layer 
analysis,  even  yielding  to  hand  computation  if  desired. 
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I.  INTRODUCTION 


It  is  the  purpose  of  this  report  to  develop  and  illustrate  a  new 
type  of  approximate  method  for  calculating  the  skin  friction  distribu¬ 
tion  in  a  compressible  turbulent  boundary  layer  under  arbitrary  heat 

transfer  and  pressure  gradient  conditions.  This  method  is  an  extension 

* 

of  an  incompres s ible  flew  analysis  reported  by  White  ( 74) . 

One  must  review  the  present  status  of  compressible  turbulent 
boundary  layer  calculation  in  order  to  justify  the  need  for  a  new 
method.  A  great  many  workers  are  active  in  the  field  of  boundary 
layer  prediction.  For  incompressible  turbulent  flow,  over  sixty 
different  methods  exist,  and  the  relative  merits  of  some  twenty- 
eight  of  these  were  thrashed  out  thoroughly  in  the  recent  Stanford 
Conference  edited  by  Kline  et  al  ( 38} .  Attention  has  now  turned  to 
the  compressible  turbulent  boundary  layer,  and  a  review  last  year  by 
Beckwith  ( 4  )  of  the  new  methods  in  this  field  contains  one  hundred 
references.  Also,  in  1968,  an  entire  symposium,  edited  by  Bertram  (5  ), 
was  devoted  to  the  compressible  turbulent  boundary  layer. 

Following  Beckwith  ( 4  ) ,  we  may  divide  the  compressible  flow 
methods  into  three  types: 

1.  finite  difference  (FD)  methods  which  attack  the  full  boundary 
layer  equations  by  dividing  the  flow  field  into  a  two-dimensional 
mesh. 

2.  Integral  (IM)  methods  which  utilize  the  compressible  form  of 


* 

Numbers  in  parentheses  denote  references,  which  are  grouped  alphabetically 
at  the  end  of  this  report. 
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the  so-called  Karman  integral  relations  -  cf.  Schlichting  (62), 
Eqs.  (13.80)  and  (13.87)  -  plus  suitable  auxiliary  information 
about  the  behavior  of  the  various  integral  thicknesses  and 
shape  factors. 

3.  Correlation  techniques  (CT)  which  relate  skin  friction  and 
Stanton  number  to  local  flow  parameters  through  empirical 
algebraic  expressions,  mostly  derived  from  flat  plate  data 
but  often  employed  in  more  general  problems. 

To  these  we  may  add  a  fourth  type  of  method  to  which,  presumably,  the 
present  report  belongs: 

4.  Methods  which  utilize  a  markedly  different  point  of  view 
without  sacrificing  either  physical  realism  or  computational 
accuracy. 

All  of  these  methods  of  course  suffer  from  the  fact  that  the  turbulence 
terms  are  not  well  defined  and  certainly  not  known  exactly  in  any  situ¬ 
ation,  particularly  for  compressible  flow.  Thus  all  computations  of 
the  turbulent  boundary  layer  are  approximate  and  semi -empirical,  even 
if  the  most  sophisticated  finite  difference  techniques  are  used.  This 
is  not  to  downgrade  the  FD  methods,  which  are  time  consuming  but 
sufficiently  comprehensive  to  allow  one  to  make  "numerical  experiments" 
into  the  nature  of  the  turbulence  approximations. 

To  resolve  the  turbulence  terms  and  achieve  closure  of  the  basic 
equations,  various  approaches  can  be  taken  for  compressible  flow: 

1.  Eddy  viscosity,  eddy  conductivity,  and  turbulent  energy 
correlations  -  primarily  used  in  the  FO  methods. 
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2.  Empirical  correlations  between  skin  friction,  Stanton  number, 
integral  thicknesses,  and  shape  factors  -  primarily  for  IM 
methods. 

3.  Compressibility  transformations  which  relate  the  compressible 
flow  to  a  supposedly  "equivalent"  incompressible  flow  -  useful 
in  all  three  types  of  methods  (FD,  IM,  and  CT). 

4.  The  law-of-the-wall  or  the -wall -of -the -wake.  These  two  laws 
are  well  accepted  physically  and  useful  in  all  methods,  in¬ 
cluding  the  present  report,  which  utilizes  the  law-of-the-wall 
as  a  sort  of  "equation  of  state"  of  turbulence.  Kline  (38) 

has  stated  that  no  method  which  ignores  the  law-of-the-wall  can 
be  successful. 

We  may  list  by  type  the  following  boundary  layer  methods  which  have 
been  applied  successfully  to  compressible  flow,  at  least  for  adiabatic 
walls : 

1.  FD  methods:  Herring  and  Mellor  (30) ,  Cebeci,  Smith  and  Mosin- 
skis  (14),  Patankar  and  Spalding  (58),  Fish  and  McDonald  (26), 
Bushnell  and  Beckwith  (11),  Bradshaw  and  Ferriss  (6  ),  Cebeci  (13). 
Computer  listings  are  available  to  the  general  reader  from  Herring 
and  Mellor  and  from  Patankar  and  Spalding. 

2.  IM  methods:  Reshotko  and  Tucker  (59),  Camarata  and  McDonald  (12), 
Alber  and  Coats  (2  ),  Henry  et  al  (29),  Nielsen  and  Kuhn  (55), 

Shang  (63),  Sasman  and  Cresci  (61),  Winter,  Smith  and  Rotta  (77), 

A  complete  FORTRAN  program  for  the  method  of  Sasman  and  Cresci 
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is  given  by  McNally  (47). 

3.  CT  methods:  Van  Driest  (71),  Spalding  and  Chi  (66),  Sommer 
and  Short  ( 65) ,  Eckert  ( 24)  ,  Moore  ( 52) ,  Coles  ( 18) ,  Komar 
{  39)  ,  Wilson  (  75) ,  and  Tetervin  ( 70) . 

As  mentioned,  many  of  these  methods  use  the  compressibility  trans¬ 
formations  which  relate  the  equations  to  an  incompressible  flow.  Such 
transformations  have  a  long  history  in  laminar  flow,  but  in  turbulent 
flow  they  were  apparently  first  suggested  by  Van  Le  (73).  Subsequently, 
the  idea  was  developed  by  Mager  ( 44) ,  Baronti  and  Libby  ( 3  ) ,  Lauf er 
(41),  Lewis,  Kubota  and  Webb  (43),  and  culminating  with  an  extensive 
recent  discussion  by  Economos  (25).  While  the  compressibility  trans¬ 
formations  are  indeed  accurate  for  flat  plate  conditions  at  moderate 
Mach  numbers  and  heat  transfer  rates  (and  the  present  report  leads 
coincidentally  to  just  such  a  transformation) ,  they  are  based  upon 
a  kinematic  invariance  between  compressible  and  incompressible  tur¬ 
bulence.  Thus  the  transformations  probably  fail  at  conditions  of 
strong  pressure  gradient,  high  Mach  number,  or  large  heat  transfer. 

The  present  method  chooses  not  to  rely  upon  such  a  transformation 
except  for  flat  plate  conditions  wher*  the  idea  arises  implicitly. 

The  present  method  is  intended  to  coqpete  with  (or  even  replace) 
the  other  integral  methods  now  in  use.  Let  us  therefore  discuss  these 
other  methods.  To  the  authors'  knowledge,  all  integral  schemes  now 
in  use  for  the  compressible  turbulent  boundary  layer  have  their  roots 
in  the  Karman  integral  relation.  This  relation  arises  by  integrating 
the  two-dimensional  compressible  momentum  equation  with  respect  to  y 
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across  the  entire  boundary  layer.  One  form  of  the  result  is  as  follows: 


0  dU  U  dp 

U  3x  t2  *  H12  *  p  dU 

6  a  e 


]  + 


where  subscript  "e"  denotes  freestream  conditions.  The  momentum  thick* 

* 

ness  0  and  displacement  thickness  o  have  their  compressible  forms: 


,  p  u  u  »  . 

{  i  V  U  -  *  >  Sr 


* 


and  the  shape  factor  ■  4*/*.  Equation  (1)  is  a  rather  general  fom 
of  the  Karman  integral  relation,  as  discussed  by  H.  McDonald  in  Bertram 
(  5 ).  it,  for  example,  the  freestreaa  is  adiabatic  -  which  is  the 
usual  case  *  the  term  involving  the  freestreaa  density  variation  nay 
be  rewritten  as: 
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(3) 


Also,  the  third  term  on  the  left,  involving  the  pressure  variation 
across  the  boundary  layer,  is  neglected  in  aost  integral  analyses. 

For  incompressible  flow,  it  is  certainly  true  that  p  *  p#  to  good 
approximation,  and  this  term  vanishes.  However,  for  compressible  flow  * 
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particularly  with  large  streamwise  pressure  gradients,  this  term  may 

be  quite  significant.  Michel  (49)  found  in  an  e^qperiment  at  M  =  2.0 

6 

that  the  wall  pressure  could  be  as  much  as  25%  higher  than  p  and 

e 

that  Eqn.  (1)  could  not  be  balanced  to  within  30%  of  the  measured 
momentum  thickness  without  the  inclusion  of  the  pressure  variation 
term,  similar  results  are  reported  by  Hoydysh  and  Zakkay  ( 35) .  It 
appears  that  integral  methods  which  neglect  this  effect  have  simply 
delayed  the  moment  of  truth  by  masking  the  error  in  a  pseudo-corre¬ 
lation  between  0,  H^2 ,  and  which  temporarily  accounts  for  the 
discrepancy.  McDonald  goes  on  to  state  that  no  Karman  integral 
method  should  neglect  the  pressure  variation  term,  and  Myring  and 
Young  (  54)  have  indicated  a  "Mach  wave"  approximation  for  calculating 
this  term. 

The  Karman  integral  relation,  then,  hardly  stands  on  its  own  as 
an  analytical  tool  for  the  compressible  turbulent  boundary  layer.  It 
is  one  equation  in  four  unknowns:  1)  0 >  2)  H^j  3)  Cfi  and  4)  the 

pressure  variation  term.  Therefore  it  must  be  liberally  supplemented 
by  other  empirical  and  analytical  relations.  The  new  relations  often 
bring  in  new  variables  and,  before  closure  is  finally  achieved,  the 
final  package  of  equations  can  be  quite  imposing.  For  example,  the 
recent  integral  method  of  Alber  and  Coates  ( 2  ) ,  which  is  one  of  the 
most  accurate  to  date,  uses  the  following  relations  1 

1.  The  Karman  integral  relation  -  Bqn.  (1). 

2.  The  mean  energy  integral  relation. 
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3.  The  law -of -the -wall. 

4.  The  law -of -the -wake . 

5.  The  generalized  velocity  variable  suggested  by  van  Driest  (71). 

6.  A  correlation  for  the  equilibrium  dissipation  integral. 

7.  An  empirical  expression  relating  wake  function  to  local  pressure 
gradient. 

8.  A  modified  Crocco  relation  for  density  variation  across  the 
layer. 

9.  The  lateral  pressure  gradient  term  in  Eqn.  (1)  is  neglected. 

Even  with  this  formidable  package  of  approximations  and  auxiliary  re la- 
tic  ns  ,  the  method  of  Alber  and  Coates  is  valid  only  for  adiabatic  flow  and 
is  yet  to  be  extended  to  heat  transfer  conditions.  Similarly,  other 
integral  methods  grow  to  substantial  slse  when  compressibility,  heat 
transfer,  and  pressure  gradient  have  all  been  accounted  for.  The 
cosputer  program  of  McNally  (47)  for  the  method  of  Sasman  and  Cresd  (61) 
contains  over  one  thousand  lines  of  FORTRAN  instructions.  This  is 

the  same  order  of  -omplexity  as  the  FD  methods,  although  in  fact  the 
integral  methods  were  intended  to  be  an  order  of  magnitude  simpler 
th*u  finite  difference  computations,  in  the  present  state  of  integral 
methods,  then,  the  working  relations  are  too  complicated  to  allow  for 
hand  computation,  and  the  canned  programs  offered  to  the  user  are 
too  explicated  for  trouble  shooting.  The  user  is  left  even  more 
impotent  by  the  FD  methods,  which  must  be  accepted  at  their  face  values 
the  product  of  years  of  work  by  their  authors.  TWo  years  ago,  the  pre¬ 
sent  writers  obtained  a  FORTRAN  deck  (800  cards)  of  one  of  the  better 
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finite  difference  methods.  The  program  runs  beautifully  with  the 
sample  data  included  in  the  instructions,  but  numerical  overflow  always 
occurs  when  the  writers'  new  data  is  inserted.  Mo  doubt  the  writers 
are  at  fault,  but  unresolved  human  failings  are  often  the  result  when 
large  computer  listings  are  borrowed  and  put  to  new  use. 

It  is  the  purpose  of  this  report  to  present  an  alternate  inte¬ 
gral  method  which  is  quite  frankly  meant  to  compete  with  the  Karman 
integral  approach.  This  may  well  be  an  impossible  task.  The  Kar¬ 
man  integral  relation  is  an  absolute  monarch  at  the  present  time.  Some 
idea  of  its  pervasiveness  throughout  the  field  of  boundary  layer 

phenomena  can  be  had  by  studying  the  recent  data  of  Brott  et  al  ( 9 ) 

* 

for  a  favorable  pressure  gradient  at  about  *  4.  After  analysing 
all  of  their  data,  these  authors  suggest  the  following  empirical  formula 
for  calculating  the  skin  friction  coefficient  C{  in  a  supersonic  flow 
with  favorable  pressure  gradient! 

Cf  *  a(B)  Reb(p)  ,  where  A  »  -  *  &  and  %  •  V/ve*  (4) 

Mere  a  and  b  are  curve-fit  functions  of  the  pressure  gradient  parameter  A, 

which  is  a  variation  of  Clauser's  (17)  original  parameter  that  wae 
* 

based  upon  o  instead  of  9.  Now  this  formula  is  dimensionally 
impeccable  and  agrees  well  with  Brott 's  measured  skin  friction,  but 
is  totally  frustrating  to  the  present  writers  and  has  little  more 
than  nuisance  value.  The  reason  is  that  both  of  the  chosen  parameters 


* Brott 'r  data  will  be  compared  in  this  report  with  the  present  method. 
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are  proportional  to  the  local  momentum  thickness  0(x) ,  which  is  un¬ 
known  apriori.  If  one  knew  0(x) ,  it  would  appear  that  one  has  actually 
solved  the  problem,  so  that  Eqn.  (4)  would  not  even  be  needed.  This 
dilemma  vanishes  if  one  adopts  an  orthodox  stance,  in  which  case  he 
is  expected  to  compute  9<x)  from  Eqn.  (1)  and  its  auxiliary  relations. 
Thus  Eqn.  (4) ,  which  is  accurate  and  concise  and  tempts  one  to  file 
it  away  for  immediate  use,  actually  has  no  intrinsic  value:  it  is 
merely  another  auxiliary  relation  for  the  Karman  integral  equation. 

In  the  present  view,  it  is  a  frustrating  and  ever  recurring  pattern 
of  correlating  Cf (x)  with  an  integral  thickness  such  as  0(x)  and  hence 
replacing  one  unknown  by  another.  The  method  to  be  presented  here 
attacks  Cf(x)  directly  and  ignores  integral  thicknesses  and  shape 
factors,  which  can  be  calculated  later  by  algebraic  formulas  if  one 
so  desires. 

The  new  method  will  be  shown  to  be  reasonably  accurate  when  com¬ 
pared  with  data.  Indeed,  it  may  be  the  most  accurate  overall  of  any 
method  yet  proposed  for  the  compressible  turbulent  boundary  layer. 

It  makes  uss  only  of  the  equation*  of  motion  plus  a  single  extra 
relation  needed  for  closure:  the  law -of -the -wall.  The  law-of -the -wake 
in  its  standard  form  is  not  used  explicitly  but  is  implied  to  the 
extent  that  deviations  from  the  logarithmic  law-of -the -wall  may  be 
called  a  "wake*.  Apparently  the  first  serious  attempt  to  develop 
this  new  idea  ms  a  paper  by  Brand  and  Person  (7)  in  1964,  con corning 
incompressible  flow  at  very  modest  pressure  gradients.  Later  work  for 
strong  pressure  gradients  was  reported  by  White  (74) .  The  present  re¬ 
port  is  the  first  attempt  to  extend  the  method  to  turbulent  compressible 
flow. 
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II.  DEVELOPMENT  OF  THE  NEW  METHOD 


The  present  analysis  is  restricted  to  steady  two-dimensional  flow 
of  a  perfect  gas  in  a  compressible  turbulent  boundary  layer.  These 
restrictions  are  not  critical,  and  the  success  of  the  method  in 
arbitrary  applications  will  govern  whether  further  generalization  is 
warranted.  Thus  we  are  concerned  with  the  following  four  basic  re¬ 
lations  for  the  turbulent  boundary  layer  mean  flow: 

a)  The  continuity  equation: 

~<0U)  ♦  *  0  (F.) 

b)  The  momentum  equation: 


o  u 


)u 
3  x 


0  v 


H 

*7 


(6' 


c)  The  energy  equation: 

0  u  ♦  o  v||«>  9  q  +  u  t  )  (7) 

d)  The  perfect  gas  lew: 

P  ■  0*T  ,  or:  T/Tw  •  ojo  .  (8) 

2 

Here  hQ  ■  h  ♦  u  /2  is  the  stagnation  enthalpy,  and  the  symbols  q  and  T 
represent  thx  heat  flux  and  shear  stress,  respectively: 

t  *  u\y  ”  0  q  •  k||*  *  o  k1*’  <93 
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Note  that  these  equations  neglect  the  lateral  pressure  gradient  dis¬ 
cussed  earlier.  Although  the  lateral  gradient  strongly  influences 
the  momentum  thickness,  as  mentioned,  it  appears  to  have  a  negligible 
effect  on  the  wall  shear  stress  approach  used  here.  Equations  (5) 
through  (9)  were  apparently  first  assembled  by  Young  (78)  and  are 
now  standard  to  compressible  turbulent  flow  analyses.  There  are 

six  unknowns  (  p ,u,v,h  ,q,  t)  and  only  four  equations  (5-8),  so 

o 

that  further  relations  are  needed.  The  standard  point  of  departure 
for  an  FD  method  is  to  correlate  the  variables  q  and  t  with  local 
conditions  into  two  additional  expressions  for  eddy  viscosity  and  eddy 
conductivity.  The  standard  IM  method  approach  is  to  combine  Eqns.  (5) 
and  (6)  and  integrate  with  respect  to  y  across  the  entire  boundary 
layer,  thus  obtaining  Eqn.  (1) ,  the  Kerman  integral  relation.  Addi¬ 
tional  relations  for  IN  methods  are  then  brought  in  as  correlations 
between  integral  parameters. 

The  present  method  chooses  two-  and  only  two-  additional  relations! 

e)  The  compressible  law -of -the -wall > 

u+  -  u/v*  -  fcn(  r*/\  ,  ,  compressibility)  (XO) 

.  H 

where  v*«(I)  is  the  friction  velocity  related  to  wall  density.  The 

F  * 

exact  form  of  Eqn.  (10)  will  follow  later. 

f )  Croc co* s  assumption  for  the  energy  equation  -  cf .  Schllchting  (62) t 

2 

Prandtl  number  *  unity  i  T*a+bu+cu.  (11) 

The  Crocco  assumption,  which  appears  to  be  quite  reasonable  for  arbitrary 
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12 


a 


(13) 


v 

w 

T  V* 


This  is  a  very  convenient  parameter,  being  directly  related  to  shear 
stress  without  integral  thicknesses,  and  was  suggested  by  the  work 
of  Mellor  (48). 

Then  the  fundamental  assumption  of  this  paper  is  that  the  velocity 
profiles  are  correlated  by  the  law-of -the -wall  in  terms  of  the  local 
skin  friction  and  the  parameters  for  pressure  gradient,  heat  transfer, 
and  compressibility  just  discussed: 


u(x,y)/v*(x)  -  u+  *  fcn(  y+,  a ,  g , Y )  (14) 

where  y+  =  yv*/vw- 

The  actual  functional  relationship  is  a  minor  detail  and  will  be  developed 
in  the  next  section.  Equations  (12)  and  (14)  provide  the  necessary 
closure  for  the  turbulent  boundary  layer  relations,  Eqns.  (5-8).  No 
further  relations  are  needed,  and  we  will  now  derive  a  single  ordinary 
differential  equation  for  computing  the  wall  skin  friction  Cf(x)  in  an 
arbitrary  compressible  turbulent  boundary  layer. 

We  should  note  that  Eqn.  (14)  is  only  an  approximation,  it  is 
believed  to  be  accurate  under  all  but  the  most  strenuous  conditions 
such  as  a  sudden  change  in  pressure  gradient  or  a  discontinuity  in  wall 
temperature.  Other  integral  methods  are  also  remiss  in  this  respect. 
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with  FD  methods  somewhat  better  in  their  ability  to  respond  to  sudden 
local  changes.  It  is  interesting  that  the  Karman  methods  begin  with 
an  exact  integral  relation,  Eqn.  (1),  and  then  introduce  approximations, 
whereas  the  present  method  begins  with  two  approximations,  Eqns.  (12,  14), 
and  then  proceeds  exactly  from  then  on. 

Let  us  now  combine  Eqns.  (5-8,  12,  14).  Equation  (5)  is  satisfied 
identically  by  the  compressible  stream  function  defined  such  that 


3  ip 

iy  =  p 


3*. 
i  x 


-  p  v 


(15) 


Note  that  ip  has  dimensions  of  viscosity.  Equation  (15)  implies  that 
the  dimensionless  stream  function  may  be  correlated  exactly  like  u+: 


<|//uw  =  fcn(  y+ ,  a,  8,  y) 


(16) 


The  independent  variables  are  now  changed  from  (x,y)  to  (x,y  ).  By 
combining  Eqns.  (14)  and  (15)  with  the  momentum  relation,  Eqn.  (6), 
we  obtain: 


V*  u+  —  v*u+  ) 
Q  X 


3  x 


V*  3  /  A 

~  ( v^u  ) 


Vw  3y+ 


dp  vft  3  r_ 

dx'  v  3y+ 

w 


(17) 
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The  differentiation  with  respect  to  y+  is  left  untouched,  but  the 
x -derivatives  are  carried  out  using  the  chain  rule: 


1-  =  H+ J.  3a  _1  ifi  J.  lx  J_ 

;)X  3X  3y+  3x3  a  +  3X3  0  +  3X  3  y 


(18) 


The  density  variation  p(x,y)  in  Eqn.(17)  is  related  to  u+  through 
the  Crocco  relation,  Eqn.(12),  and  the  pressure  gradient  is  related  to 
the  freestream  conditions  from  the  Bernoulli  inviscid  relation: 


&  . 

dx 


duP 

“  P  u  T- 
Me  e  dx 


(19) 


By  introducing  Eqns.(18)  and  (19)  and  the  definitions  of  (  n  ,  ,  y  ) 

into  Eqn. (17) ,  we  obtain: 


p  v* 


o  +  )U  . 

Bu  ,5-3  + 


2YU+  H  * 


3J(  3u+ 
pv^  3  8  3y+ 


3  y  3y+ 


+ 


p  v* 


2  da 
dx 


1  3J;  3u\ 

pv^  3  a  3y+' 


dUe 


U 

e  e  dx 


v* 

v 

w 


3  T 
3y+ 


(20) 


Finally,  introduce  the  density  from  Eqn. (12)  and,  following  White  (74), 
perform  the  critical  step:  integrate  the  entire  Eqn. (20)  from  y  *  0 
{  t  ■  t  )  to  y  ■  6  (  t  -0) .  The  result  is  the  following  ordinary 
differential  equation  for  the  friction  velocity: 
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G  + 


i 


(21) 


dv* 

<^v*  5 


t  v!  Ss  H 
Mw  dx 


n  Me 
pe  e  dst 


v*  T, 


y 


w 


where  is  the  value  of  y+  at  ( y=S) .  The  functions  G  and  H  are 
short  notation  for  rather  lengthy  (but  straightforward)  integrals 
involving  u+  and  the  law-of -the -wall  parameters: 


G  =  /  £(u+2-6u+lii  +2yu+lii  + -ft.il  12. -3L.|iU)  dy+ 

0  Pw  y  pv)$?v+  e 


(22) 


-  /  T  (U 


+  )U 

rs 


_JL  l±  )  dy+ 


(23) 


These  expressions  will  be  evaluated  numerically  in  the  next  section. 

After  considerable  inspection,  it  may  be  verified  that  3qn.  (21) 

contains  only  a  single  unknown:  the  wall  shear  velocity  v*.  The 

remaining  quantities  such  as  G,  H,  etc.,  can  be  directly  correlated, 

through  our  assumed  law-of -the -wall  relation,  Eqn.  (14),  to  v*  and  the 

known  freestream  and  wall  temperature  conditions.  Thus  Eqn.  (21) 

is  self-sufficient:  a  first  order,  nonlinear,  ordinary  differential 

equation  in  v*(x),  which  needs  only  the  single  initial  condition  v*  • 

v*  at  .x  ■  x  .  It  is  convenient  to  non-dimensionalise  everything.  Let 
o  o 

L  and  u  be  a  reference  length  and  velocity,  respectively,  and  define 
o 

the  following  dimensionless  variables: 
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x*  *  j“  )  V  =  U/tJ  -  V(x*)  ;  X  *  (2/C.)>*  (24) 

Leo  £ 


Then  Eqn.  (21)  may  be  rewritten  in  the  following  form: 


where  the  primes  signify  differentiation  with  respect  to  x*.  Equation 

(25)  is  the  central  result  of  this  report.  It  is  valid  as  an  approximate 

means  of  solving  for  the  skin  friction  distribution  X  (xM  for  any 

arbitrary  freestream  Mach  number  and  wall  teaperatvre  distribution. 

The  proper  boundary  c  sndition  is  a  single  known  value  X  -  X  at  x*  - 

o 

x*Q.  The  only  other  relations  needed  are  1)  suitable  formulas  for 
the  quantities  G,  H,  and  5  as  functions  of  (  «,B,y,X)»  and  2)  known 
variations  with  x*  of  the  freestream  velocity  U#  and  the  wall  temperature 
ratio  T^/T^.  ihe  effective  Reynolds  number  is  defined  ass 
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The  viscosity  ratio  in  Eqn.  (26)  can  be  evaluated  in  terns  of  (T  /r  ) 

e  w 

through  any  convenient  formula  (Sutherland  law,  etc.).  For  our 
purposes,  the  simple  power-law  expression  is  quite  accurate: 


>A  *  <W 


<271 


and,  in  the  computations  which  follow,  the  value  of  n  =  0.67  for  air 
was  used.* 

For  use  in  Eqn.  (25) ,  the  pressure  gradient  parameter  a  may  be 
written  in  terms  of  the  new  variables  as  follows: 


(1/V) * 


(28) 


and  formulas  for  G  and  M  will  be  confuted  from  Eqna.  (22)  and  (23). 

The  final  thread  in  the  fabric  is  the  evaluation  of  the  thickness 
function  6+,  which  happens  to  cause  the  only  algebraic  difficulty  in 
the  entire  analysis.  Xt  is  computed  by  evaluating  Eqn.  (14)  at  y  •  6 : 


0#/v* 


x<T#/rw) 


fcn(  6  »  o»  | »  y) 


(29) 


ft 

The  commonly  used  value  n  “  0.76  is  actually  not  very  accurate  for  air. 


Hence  6+ is  implicitly  a  function  of  X  and  known  functions  of  x*. 

The  authors  have  not  been  able  to  invert  Eqn.  (29)  explicitly  for 
6+  and  thus  have  had  to  settle  for  an  iterative  procedure  of  com¬ 
puting  6+when  X  is  known.  For  the  computer  program  listed  in  the 
appendix,  this  iteration  is  no  trouble  whatever,  but  obviously  it  would 
create  tedium  in  a  hand  computation.  Eqn.  (25)  is  entirely  amenable 
to  hand  confutation  or  even  graphical  analysis  (since  it  is  only  first 
order),  but  in  such  cases  a  set  of  charts,  detailing  the  velocity  re¬ 
lationship  defined  by  Eqn.  (29) ,  would  be  handy  for  computing  . 

One  remark  is  in  order:  the  basic  relation,  Eqn.  (25),  has  the 
unique  property  among  IM  methods  of  providing  an  explicit  flow  separation 
criterion,  which  occurs  when  G  *  3  aH .  This  is  so  because  that  event 
will  result  in  X  approaching  infinity  and  hence  Cf  approaches  sero, 
which  is  the  precise  definition  of  "separation".  Thus,  in  the  present 
method,  we  need  not  search  for  a  pseudo  warning  of  separation  such  as 
a  particular  value  of  the  shape  factor  or  otherwise.  This  seems 
a  distinct  advantage,  since  the  writers  know  of  no  other  effective 
separation  criterion  for  compressible  flow.  Finally,  we  may  note  that 
Eqn.  (25)  is  identical  in  form  to  the  incompressible  analysis  of  white  (74). 
As  8  and  y  become  very  small  (negligible  heat  transfer  and  compressi¬ 
bility)  ,  Bqns.  (22)  and  (23)  for  G  and  R  become  identical  to  that 
earlier  analysis. 
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III.  THE  COMPRESSIBLE  FLOW  LAW -OF -THE -WALL 


To  complete  the  analysis  and  make  Eqn.  (25)  useful,  we  must  develop 
a  quantitative  formula  for  the  law-cf -the -wall  as  defined  by  Eqn.  (14). 

The  authors  tried  many  approaches,  the  mo3t  promising  of  which  were? 

1.  the  effective  velocity  concept  of  van  Driest  (71); 

2.  the  transformation  theory  of  Coles  (18); 

3 .  che  transformation  of  Baronti  and  Libby  ( 3 ) ; 

4.  correlation  of  measured  profiles,  e.g.  Kepler  and  O'Brien  (37); 

5.  the  eddy  viscosity  approach  of  Deissler  (23). 

The  van  Driest  effective  velocity  was  used  in  the  integral  method  of 
Alber  arid  Coats  (  2  ).  Maise  and  McDonald  (45)  found  that  adiabatic 
supersonic  boundary  layer  profiles  would  collapse  fairly  well  (+30%)  to 
the  incompressible  law  of  the  wall  u  (y  )  if  the  local  velocity  u  were 
replaced  by  the  van  Driest  (71)  generalised  velocity  u*t 

u*  »  (U  /alsin^fau/U  )  ,  where  a2  »  rM  /(1+rM  ).  (30) 

e  a  e  • 

The  concept  can  be  extended  to  flow  with  heat  transfer,  also,  but  the 
agreement  is  much  poorer,  especially  for  high  Mach  numbers. 

The  various  transformation  theories  -  see  Sconomoe  ( 25)  for  a  recent 
review  -  are  becoming  lass  popular  because  of  their  algebraic  complexity 
and  because  they  stand  on  very  weak  ground  in  pressure  gradient  or  heat 
transfer  conditions.  Nevertheless,  the  approach  is  viable,  and  the  trans¬ 
formations  of  Coles  ( 18)  and  of  Baronti  and  Libby  ( 3  )  have  recently  been 


suggested  by  Hopkins  et  al  (34)  as  a  means  of  correlating  flat  plate 

skin  friction.  However,  every  single  computation  of  a  given  point 

■f  4 

u  (y  (a  j  y)  by  transformation  theory  requires  multiple  iterations, 
and  for  the  present  application  the  idea  was  finally  dropped. 

Scheme  number  four,  the  correlation  of  measured  velocity  profiles, 
is  an  obvious  resort  for  finding,  at  least  roughly,  the  effect  of  the 
parameters  (a,6,Y)  on  the  law-of-the^wall.  Kepler  and  O'Brien  (37) 
measured  supersonic  adverse  pressure  gradient  profiles  (positive  a  and 
Y  ),  while  Lee  et  al  (42)  measured  cold  wall  heat  transfer  (positive  8  ) 
and  Brott  et  al  ( 9 )  studied  supersonic  favorable  pressure  gradients 
(negative  0  ,  positive  Y  ).  From  these  we  find  that  the  general 
trend  of  effects  is  as  sketched  in  Figure  1.  The  effect  of  positive 

<*  and  positive  8  is  to  raise  the  date  above  the  familiar  incompressible 
logarithmic  lew,  and  negative  a  and  8  have  the  opposite  effect.  The 
parameter  Y  ,  being  proportional  to  velocity  squared,  is  always  positive 
and  always  tends  to  lower  the  data  as  shown.  If  (o,8,y)  are  all  finite, 
the  general  effect  is  roughly  a  superposition  of  the  various  separate 
effects.  There  is  such  a  large  scatter  in  measured  supersonic  profiles 
and  skin  friction  that  no  quantitative  effects  could  be  correlated. 

The  final  scheme  is  an  eddy  viscosity  approach,  using,  say,  a 
mixing  length  approximation  and  the  Crocco  relation  for  temperature. 

This  approach,  which  is  simple  and  fairly  quantitative  and  allows 
ready  evaluation  of  the  integrals  in  Kqns.  (22)  and  (23),  was  adopted 
for  the  present  study.  It  agrees  with  the  effects  sketched  in  Figure  1 
and  gives  numerical  values  of  the  functions  G,H,  and  6twhich  seen  quite 
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Figure  1.  SKETCH  ILLUSTRATING  THE  EFFECT  OF 

VARIOUS  PARAMETERS  ON  THE  LAW- OF- THE- WALL 
FOR  COMPRESSIBLE  TURBULENT  FLOW. 


adequate  for  general  use,  as  shown  in  the  next  two  sections. 

The  theory  chosen  was  that  of  Deissler  (  23) ,  who  assumed  that 
the  Prandtl  mixing  length  approximation  could  be  extended  to  the  variable 
density  case  with  no  further  changes.  That  is,  the  total  shear  is  re¬ 
lated  to  an  eddy  viscosity  as  follows s 


t  »  (u  +  c)  fir  ' 

**  (31) 


.  2  2  |  »u  | 
where  e  -  p  ^  y  |  ly  |  ' 

and  k  *  0.4  is  von  Karman’s  constant.  Also,  near  the  wall,  the  boundary 
layer  is  approximately  a  Couette  flow  with  negligible  convective  accelera¬ 
tions,  so  that  the  egression 


*  4  t +  f?  y  *  or*  r  *  i  +  «  y+  (32) 

w  dx  Tw 

is  a  good  approximation  and  used  in  many  theories.  Finally,  the  density 
in  Eqn.  (31)  is  eliminated  through  the  perfect  gas  relation,  Eqn.  (8). 
Equations  (31)  and  (32)  may  be  combined  to  yield 

i*  (  *2  ***>  >£'.  >  '  <»> 

which  is  quadratic  in  the  velocity  derivative.  Solving,  with  viscosity 
and  density  replaced  by  tmperature  through  Kqns.  (UJ  and  (27),  we  obtain t 
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(34) 
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Equation  (35)  may  then  be  integrated  -  using  e.g.  Subroutine  RUNGE  in 
the  appendix  -  to  obtain  numerical  values  of  u+(y+,a»8»y) •  Closed 
forms  for  the  integral  can  also  be  found,  but  the  resulting  expressions 
are  cumbersome  algebraically.  He  may  note  that  Eqn.  (35)  is  precisely 
the  expression  used  in  the  finite  difference  method  of  Patankar  and 
Spalding  (58)  for  confuting  local  velocity  profiles  near  the  wall. 

Some  numerical  values  of  the  integral  of  the  Eqn.  (35)  are  shown  in 
Figures  2  and  3.  For  (a,(S,y)  =  (0,0,0),  the  result  is  the  familiar 
logarithmic  l*.w: 

u+  -  ^  ln(y+)  +  5.5  (36) 

For  finite  (  the  variables  can  be  separated  in  Eqn.  (35)  and 

integrated  to  give  the  following: 

(37) 

2  >1  + 

where  Q  =  (8  +  4  y)  and  P  =  (1  +  a  Jr  )  • 

Since  this  formula  diverges  at  (y  =  0) ,  it  is  necessary  to  insert 
initial  values  (u*,  y*)  sufficiently  close  to  the  wall  for  all  the 
curves  to  converge.  From  Figures  2  and  3,  this  appears  to  happen  for 
y+  less  than  about  ten.  Thus,  for  utilizing  Eqn.  (37)  in  calculations, 
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WALL  FROM  EQ.  (37)  FOR  ZERO  HEAT  TRANSFER 


Figure  3.  THE  LAW-OF-THE- WALL  FROM  EQ.  (37)  FOR  FINITE  HEAT  TRANSFER. 


+  + 

we  will  take  as  initial  conditions  the  point  (uq  =  10,  yQ  ■  6),  which 
falls  on  the  logarithmic  curve,  Eqn.  (36).  Note  that,  if  y+  is  known 
u+  may  be  computed  explicitly  from  Eqn.  (37) ,  but  the  opposite  is  not 
true,  as  mentioned  earlier.  Unfortunately,  it  is  the  problem  of  com¬ 
puting  y+  when  u+  is  known  that  confronts  us  in  the  present  analysis 
(for  evaluating  6+in  Eqn.  (25)),  hence  the  need  for  iteration. 

Figure  2  shows  the  law-of -the -wall  profiles  for  a  range  of  values 
of  a  and  y  ,  while  Figure  3  shows  the  effect  of  f?  .  The  trend  is 
exactly  as  illustrated  in  Figure  1.  No  attempt  will  be  made  to  improve 
upon  these  curves  by,  say,  including  a  "wake".  As  sketched  in  Figure  4, 
Eqn.  (37)  will  simply  be  assumed  to  hold  all  the  way  to  the  edge  of 
the  boundary  layer  whereas  the  dotted  line  illustrates  the  true  outer 
wake.  Thus  we  make  some  error  in  estimating  the  boundary  layer  thick¬ 
ness,  but  the  effect  on  skin  friction  is  entirely  negligible,  particularly 
for  a  supersonic  turbulent  boundary  layer,  which  shows  almost  no  wake. 
There  is  also  only  about  a  three  per  cent  error  in  using  this  approxi¬ 
mation  to  calculate  the  momentum  and  displacement  thicknesses,  if  one 
should  care  to  compute  such  quantities.  Note  in  passing  that  the 
arcsines  in  Eqn.  (37)  can  never  have  an  invalid  argument,  because  of 
the  physical  requirement  that  (T/T^)  from  Eqn.  (12/  remain  positive. 
However,  if  a  is  negative  (favorable  pressure  gradient),  there  is  a 
mathematical  possibility  that  the  quantity  (1  +  oyt  could  become  nega¬ 
tive  at  large  y+,  thus  rendering  P  imaginary  and  the  formula  invalid, 
when  this  happens,  though,  the  velocity  gradient  is  negative,  implying 
that  the  velocity  inside  the  boundary  layer  has  become  greater  than 
the  frees tream  value  Uq.  To  avoid  this  difficulty,  it  is  recoimnended 
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Figure  4.  COMPARISON  OF  THE  PRESENT  NEAR-WALL 
THEORY  WITH  THE  ACTUAL  LAW-OF-THE-WAKE. 


A  third  function  of  interest,  to  be  called  F,  is  the  grouping  which 
occurs  in  our  basic  relation,  Eqn.(25): 

F(u+,  a  ,  8  ,  Y  )  =  X  (  X2  «+  -  G)  (40) 

This  quantity  also  has  a  limiting  exponential  approximation: 

F(u+, 0,0,0)  =  47.0  exp(0.475  u+)  ■  5.53  G(u+, 0,0,0)  (41) 


It  was  this  near-proportionality  of  F  with  G  that  enabled  White  (74) 
to  find  an  exact  solution  to  Eqn.(25)  for  modest  pressure  gradient 
in  incompressible  flow.  The  same  fact  will  also  enable  us  to  simplify 
our  calculations  for  the  present  study. 

By  constructing  curve-fit  expressions  which  reduce  to  the  above 
limiting  formulas  for  (o,8,y)  *  (0,0,0),  we  have  adopted  the  following 
approximations  which  appear  to  give  good  accuracy  for  the  practical 
range  of  values  of  all  three  parameters: 

G(u+,o  ,8  ,v  )  *  8,5  exp  f  0.475  u+f/ft+0.1  sgn(a)  («8+),|l) 

e  9 

H(u*,a  ,8  ,y  )  *  0.062  exp  f  0.84  u*f/ri+0.l2  sgn(«)  (««+)0*6 1) 

(42) 

F(u*,0,8,Y)  *  5.53  G  , 

e 

where  f  -  (1  ♦  0.22Y  u^J/U  ♦  0.38  u*) 
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These  approximations  are  utilised  with  Eqn.  (25)  to  compute  the  skin 

2 

friction  distribution  Cf  *  2/  x  .  The  factor  f  above  will  play  a 
prominent  role  in  the  flat  plate  theory  of  the  next  section.  The 
parameter  a  is  computed  from  Eqn.  (28)  and  u+  is  related  to  the 
variable  X  from  Eqn.  (29) .  For  a  perfect  qas  of  specific  heat  ratio 
k,  the  two  quantities  in  the  factor  f  in  Eqn.  (42)  may  be  written  as 
follows: 


Y  uj2  -  r  (^)  t?  (Te/Tw) 

8ue+  -  (Te/rw)  (1+r^ji  M^)  -  1  *  <T*/V  -  1 


(43) 


Finally,  since  the  correlations  in  Eqn.  (42)  contain  5+,  Eqn.  (37) 
must  be  iterated  for  a  given  value  of  u*  to  compute  this  thickness 
parameter.  As  mentioned  earlier,  the  computation  of  6+  is  the  only 
cumbersome  procedure  associated  with  the  present  method.  It  is  hope'* 
that  future  work  will  enable  6+  to  be  alimiru*  ~ed  entirely  in  favor  of  the 
single  variable  X  . 

This  new  theory  will  now  be  illustrated  in  the  next  two  sections 
for  cases  of  practical  interest.  The  complete  set  of  basic  equations 
is  listed  as  a  group  in  the  last  section  of  this  report. 


33 


IV.  COMPRESSIBLE  TURBULENT  FLOW  PAST  A  FIAT  PLATE 


Needless  to  say,  the  flat  plate  is  not  the  reason  for  developing  the  new 
theory,  even  if  heat  transfer  is  included.  There  have  been  numerous 
theories  of  compressible  turbulent  flow  over  a  Elat  plate,  and 
twenty-four  of  these  were  discussed  in  detail  in  1963  in  a  review  paper 
by  Spalding  and  Chi  (66),  who  developed  a  least-squares  correlation  of 
their  own.  since  then,  another  dozen  methods  have  appeared,  mostly  con¬ 
cerned  with  the  recent  interest  in  generalised  compressibility  trans¬ 
formations.  To  add  on  yet  another  method,  Eqn.  (25) ,  is  not  an  rtr  ng 

prospect.  Yet  a  comparison  should  be  made  for  completeness,  and  to  this 
end  the  Spalding  and  Chi  (66)  data  compilation  was  brought  up  to  date. 

The  Appendix  lists  data  for  flat  plate  flow  at  426  adiabatic  wall  condi¬ 
tions  and  147  heat  transfer  conditions,  at  freeatream  Mach  numbers  from 
zero  to  ten.  It  is  apparently  the  largest  list  ever  compiled  of  such 
data.  The  measurements  were  compared  with  the  various  theories  (see 
Table  1)  with  two  very  interesting  results:  1)  for  an  adiabatic  wall, 
the  present  method  has  the  smallest  mean  absolute  error  of  any  known 
theory:  and  2)  for  flow  with  heat  transfer,  an  almost  forgotten  method 
given  in  a  dissertation  by  Moore  (52)  is  the  clear  winner  -  the  present 
method  scoring  only  a  reasonably  creditable  fourth  place,  behind  Spalding 
and  chi  (66)  and  van  Driest  (72).  Of  critical  importance  is  the  curious 
fact  that,  of  all  tna  various  flat  plats  theories,  only  the  present 
method  car.  be  extended  without  any  complication  to  variable  frees '.ream 
and  wall  temperature  conditions.  Thus  che  writers  believe  that  they  have 
put  forth  a  demonstrably  accurate  *nd  generally  useful  new  method. 


For  the  flat  plate,  V  8  1.0  8  constant,  and  Eqn. (25)  reduces  to: 


G  d  X  8  R^  dx*  (44) 

If  the  freestream  Mach  number  and  wall  temperature  are  constant,  8  and 
Y  are  constant,  and  G  is  therefore  a  function  of  X  only.  He  may 
integrate  Eqn. (44)  to  obtain  a  relation  between  local  Reynolds  number 
and  local  skin  friction: 

X 

/  G(X.8,y)  dX  (45) 

0 

where  we  have  assumed  fully  turbulent  flow  with  the  leading  edge  (x  8  0) 
at  the  beginning  of  the  boundary  layer  (  X8  0).  The  integration  could 
be  performed  numerically,  but  we  desire  an  explicit,  albeit  approximate, 
skit,  friction  foimula.  Thus  we  adopt  the  curve-fit  expression  for 
from  Eqn.  (42),  for  <*  ■  0: 

G<  M,1»  *  8.5  exp(  0.475  f  X  <T  A  )S 

substituting  into  Eqn. (44)  and  Integrating,  we  obtain: 

0.475  f  \  <T  ,T  )**  *  ln(  0.056  f  (T  /T J1+nR  ) 

*  w  e  w  x 

15 


(46) 


(47) 


R*  8  oex/v48  (T/re> 


n+*» 


Solving  for  =2 /  \2,  we  obtain  the  following  interesting  formula: 


C 


f 


0.451  f2  (TeAw) 

In2  (0.056  f.  (T  A1  )1+nR  ) 
e  w  x 


(48) 


The  factor  f  is  defined  in  Eqn.(42)  and,  for  a  perfect  gas,  becomes: 


(  1  +  0.22  r  (^i)  M2  T  A  ) 

f  =  - 2 - i-JLJL.  (49) 

(1+0.3  (TawAw-  > 


In  the  limit  of  incompressible  flow  with  zero  heat  transfer,  f  =  1.0  and 

T  =  T  =  T  ,  and  we  obtain: 
e  w  aw 


Cf  (Incompressible  Flat  Plate)  -  — - - _  (50) 

In  (0.056  R  ) 
x 


This  formula  was  obtained  by  white  (74)  in  his  incompressible  analysis. 

Now  neither  Eqn. (48)  nor  (50)  is  particularly  accurate,  being  based  upon 
a  curve-fit  for  G,  but,  by  comparing  the  two,  we  may  deduce  the  following 
transformation  for  flat  plate  flow: 


C,(R  ,M  ,  T  A  )  = 

f  x  e  e  w 


f  -e  C.  (R  .  .  ) 

T  f.  effective 
w  me 


(51) 


where  R 


effective 


f  (T  A  )1+n  R  • 
e  w  x 


36 


Thus  the  flat  plate  skin  friction  for  compressible  flow  is  directly 

related  to  an  incompressible  value  of  Cf  evaluated  at  a  different 

(usually  smaller)  Reynolds  number.  This  is  the  same  concept  achieved 

by  previous  theories  of  the  flat  plate  in  turbulent  flow.  In  the  notation 

of  Spalding  and  Chi  (66),  C_  should  be  related  to  by  the  relation: 

rinc 


cf  =  p  cf  <RX  fr  > 

c  inc  x 


F  %  (R©  V  ' 

c  inc  © 


(52) 


depending  upon  whether  the  appropriate  known  Reynolds  number  is  or 

R  .  The  uniqueness  of  this  transformation  is  satisfied  only  if 

w 


(53) 


From  Eqn. (52),  we  see  that  the  present  theory  predicts  that 


=  !w  f-2 
Te 


f  (T  /T  )1+n  ;  F„  =  (T  /T  )“  f"1 
e  w  RQ  e  w 


(54) 


with  f  given  by  Eqn. (49).  There  are  many  other  theories,  and  Table  1 
summarises  eight  methods  which  are  either  well  known  or  very  accurate  or 
both.  Note  the  diversity  of  expressions  for  and  Fr^  .  The  first 

six  of  these  were  tested  against  the  573  data  points  in  the  Appendix.,  and 
the  relative  deviations  e^  were  calculated  according  to  the  formula: 
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TABLE  1 

FLAT  PLATE  TRANSFORMATION  FUNCTIONS 


AUTHOR 

F 

c 

NOTATION 

Eckert (24) 

T_A 

E  e 

<W" 

te=  Te(l+.039M2-.5(l-l/t>) 

t  =  T  /T 
e  w 

Moore  (52) 

(1/t-l) 

Qfsin^d-t)**)2 

.1156Z/L 

t  =  T  A 
g  w 

Q  =  .9212  e-0706^) 

L  -  11.5  +  6.6(1-T  A  ) 
w  aw 

Z  =  tn  exp(.4L) 

Sommer  & 
Short  (65) 

TSAe 

<w” 

Te=T  (1+.035M2+. 45 (1/t-l)) 

S  ©  6 

Spalding  & 
Chi  (66) 

.2  r  M2 
-1  e_  1  0 
(sir.  A+sin  B) 

.702  .772 
t  J 

t  =  T  A 
e  w 

a  =  (.2r  M2  t)*1 

b  =  t(l+.2r  M2  -  1/t) 

2  2  2  h 

A=  (2a  -b) /(4a  +b  ) 

2  2  >i 

B=  b/ (4a  +b  ) 

J  *  T  A 

w  aw 

van  Driest 
#2  (72) 

same  as 
Spalding 
and  Chi 

<T  /T  >n 
e  w 

same  as  Spalding 
and  Chi 
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TABLE  1  (Continued) 


AUTHOR 


Present 

method 


F 

c 

FRq 

NOTATION 

-1  f-2 

tn/f 

t  =  T  A 
e  w 

l+.044rM^  t 

f  1+. 3  (T  A  -1) 
aw  w 

Baronti  & 
Libby  ( 3 ) 


j  =  T  A 
J  we 


(TfAe)  (j  +  (l+.2Me 
-j)(-  3)  (^f  ±)  ^ 

-.lM2(37.45Cf .))_1 
e  l 


T  =  T  ( j+(l+.2M  -j) (10.6) 
re  e 

(Cf  (10. 6)2) 

(ITERATE  FOR  Tf ,  Cfi) 


j  =  T  A 
w  e 


To  + 

T  =  — ®  (T A  )  dy 

C  430  1  'V  y 


TA  =  j  +  (1+.2M 
e  e 

-j)  (u+)  (Cfi/2)lj 
-  u+2(.l  M2  Cfi) 


Note:  u  given  as  a 
function  of  y+  in 
reference  (18). 

ITERATE  FOR  T^,  C^. 
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Error:  e.  =  (  Cf  /  Cf  -  1)^  (55) 

exp  theory 

Both  the  rms  and  mean  absolute  error  were  computed  for  the  first  six 
methods  in  Table  1: 

1  2  **  1 

RMS  Error  =  (  -  Z  e.)  ;  Mean  Absolute  Error  =  -  £  |e.  I  (56) 

N  1  NX' 

The  final  two  methods  in  Table  1  -  the  transformation  theories  of  Coles (18) 

and  Baronti  and  Libby  (3) -were  not  computed.  A  glance  at  the  Table  shows 

why.  Not  only  do  the  compressibility  transformations  require  iteration 

to  compute  the  "reference"  temperatures  (T  ,T„,T  ),  but  within  this 

B  f  C 

iteration  is  another  nested  iteration  for  the  incompressible  skin 
friction  Cf^.  The  writers  were,  frankly,  not  prepared  for  such  an 
undertaking  and,  frankly,  doubt  if  the  average  engineer  is  prepared  for 
it  either.  Nevertheless,  the  two  theories  continue  to  be  popular. 

The  data  are  in  the  form  of  skin  friction  or  drag  coefficient  CQ 

for  various  R^  or  RQ.  Since  Fr^  and  Fr^  are  almost  always  fractions, 
the  incompressible  value  Cf ^  is  invariably  evaluated  in  the  low 
Reynolds  number  range  and  it  is  important  to  have  a  good  formula  for 
computing  Cfi  from  Reffective*  Otherwise  one  will  add  on  the  error  in 
the  formula  to  the  error  in  the  transformation.  Apparently  the  best 
incompressible  formulas  are  those  of  Spalding  and  Chi  (66),  who  integrated 
the  exact  law-of-the-wall  for  zero  pressure  gradient,  including  the 
sublayer  region.  This  results  in  implicit  formulas  for  C^: 


40 


Rx  =  75072  +  7^g(  eY(Y2-4Y+6)-6-2Y-Y4/12-Y5/20-Y6/60-y7/252)  (a) 

(57) 

2 

Re  =  756  +  O’  (  eY(1"2/Y)  +  2/y  “y2/6  “y3/12  -y4/40  “  y5/180)  (b) 

where  Y  =  0.4  ( 2/Cf . 

If  the  Reynolds  number  is  known,  one  must  iterate  these  formulas  to  find 
Cf.  by  taking  as  a  first  estimate,  say,  a  simple  power-law  formula. 
Equations  (57),  although  cumbersome,  have  the  best  known  agreement  with 
incompressible  friction  data  -  see  Spalding  and  Chi  (66)  -  and  hence  were 
used  for  the  theoretical  comparison  in  Table  2. 

If  we  wish  to  compute  the  drag  CD  for  a  given  R^,  Equation  (57a) 
is  used  to  compute  Cf(L),  which  specifies  Y (L)  and  therefore  specifies 
Rg(L)  from  Eqn. (57b).  Then  the  drag  follows  from  the  Karman  integral 
formula  for  zero  pressure  gradient: 

CD  =  2  0(L)/L  =  2  R  (L)/  .  (58) 

Again  the  computation  is  obviously  cumbersome,  but  the  accuracy  is 
desirable  for  the  theory  comparison. 

Table  2  shows  the  comparison  of  the  first  theories  from  Table  1  with 
all  of  the  flat  plate  friction  data  known  to  the  writers.  For  an  adiabatic 
wall,  the  methods  of  Spalding  and  Chi  and  of  van  Driest  are  essentially 
identical  and,  together  with  the  present  theory,  are  the  most  accurate 
theories  available.  The  present  theory  has  a  significantly  lower  mean 
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TABLE  2 


COMPARISON  OF  SIX  THEORIES  WITH  FLAT  PLATE  FRICTION  DATA 


ADIABATIC:  426  Points 

COLD  WALL:  147  Points 

AUTHOR 

RMS  %  Error 

ABS  %  Error 

RMS  %  Error 

ABS  %  Error 

Eckert  (24) 

12.44 

9.06 

24.61 

20.28 

Moore  ( 52 ) 

. 

8.87 

6.54 

12.56 

10.17 

Sommer  and 
Short  (65) 

9.40 

7.77 

19.95 

16.45 

Spalding  & 

Chi  (66) 

7.59 

5.46 

13.23 

11.04 

Van  Driest 

n  (72) 

7.55 

5.46 

14.28 

11.73 

Present 

Theory , 

Eqn. (54) 

7.59 

5.17 

14.80 

12.29 

Notes 


RMS  Error  = 


H 

) 


MEAN  ABS  Error  = 


1 

N 


E 
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absolute  error. 


For  flow  with  heat  transfer  (cold  wall  data),  a  surprising  winner 

emerges:  a  dissertation  by  Moore  (52)  which  received  only  limited 

distribution  and  was  not  included  by  Spalding  and  Chi  (66)  in  their 

review.  The  method  of  Spalding  and  Chi  takes  second  place  over  van 

772 

Driest  by  virtue  of  the  term  (T  /T  )*  added  to  Fr^.  The  present 

w  a.w  w 

theory  is  a  close  fourth  behind  that  of  van  Driest.  For  heat  transfer 
work,  the  present  theory  would  be  greatly  enhanced  by  dropping  the 
Crocco  assumption,  Eqn.(12){  its  performance  is  creditable  in  any  case. 
Note  that  the  widely  used  "reference  temperature"  method  of  Eckert  is 
in  fact  a  very  poor  performer. 

It  is  interesting  that  the  method  of  Moore  (52)  has  surpassed  the 
least  squares  data  correlation  of  Spalding  and  Chi  (66)  by  simply  coming 
up  with  a  better  formula  than  the  one  which  Spalding  and  Chi  minimized. 
Moore  used  the  van  Driest  effective  velocity,  Eqn.(30),  and  added  a 
factor  "L"  (see  Table  1)  to  account  for  the  variation  in  viscous  sublayer 
thickness  with  wall  temperature.  Also,  Moore  introduced  a  "Q"  factor 
to  modify  the  momentum  analysis  of  Wilson  (75)  to  account  for  an  adiabatic 
recovery  factor  of  r  =  0.89,  which  was  the  value  used  in  all  calculations 
for  preparation  of  Table  2.  it  appears  from  the  writers'  study  that  the 
increased  effectiveness  of  Moore's  method  is  chiefly  due  to  the  sublayer 
correction  factor  "L". 

For  the  reader's  interest,  Figure  6  shows  the  classic  (Cf/Cfi)  plot 

7 

versus  Mach  number  for  a  10  and  various  wall  temperature  ratios.  As 
is  well  known,  there  is  a  large  effect  of  wall  temperature  over  the  entire 
Mach  number  range.  Similarly,  Figure  7  shows  the  effect  of  Reynolds 
number  on  the  adiabatic  wall  friction  factor.  We  see  that  this  effect 
also  can  be  substantial  at  very  high  Mach  numbers. 
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Figure  7.  EFFECT  OF  REYNOLDS  NUMBER  ON  THE 

RATIO  OF  COMPRESSIBLE  TO  INCOMPRESSIBLE 
SKIN  FRICTION:  PRESENT  THEORY,  EQNS.  (54, 57). 
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Although  they  were  adopted  for  the  rigorous  comparison  in  Table  2, 

the  implicit  formulas  for  and  from  Eqns. (57,58)  are  really  too 

cumbersome  for  routine  use.  Therefore,  for  general  interest,  the  writers 

compared  Eqns, (57)  and  (58)  with  some  of  the  more  popular  engineering 

formulas  for  drag  and  skin  friction,  computing  the  per  cent  deviation 

from  Eqns. (57)  and  (58).  For  evaluating  (R^) ,  we  may  cite  the 

5  9 

following  formulas  and  their  accuracy  over  the  range  R^  =  10  -10  : 

a)  Blasius  power-law: 

Cf  =  0.0592/  R®*2  +  30% 

b)  Schultz-Grunow: 

Cf  »  0.370/(log1QRx)  +  8% 

c)  Prandtl-Schlichting: 

Cf  *  1/(2  log10Rx  -  0.65)2*3  +  8  %  (59) 

d)  von  Karman: 

CfmH  «  4.15  log1Q(Rx  Cf)  +  1.70  +  7  % 

e)  Present  approximation :  Eqn. (50)  modi fled - 

Cf  "  0.42/ln2 (0.056  Rx>  ♦  4) 

It  seems  that  the  present  theory  gives  the  most  reliable  approximation 
for  computing  skin  friction  explicitly. 

ror  computing  from  Rg  ,  we  have  the  following  formulast 
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a)  Prandtl-Falkner  power  law: 

Cf  «  0.013/  Re1/6  +  6  % 

b)  Squire -Young: 

Cf  *  0.0576/  109^(4.075  R0)  +  7  %  (50) 

c)  von  Karaan-Schoenherr: 

C_  *  1/(17. OS  R  2  +  25.11  Ra  +  6.012)  +  3  % 

Here  the  accuracy  is  much  better,  even  for  the  crude  formulas,  because 
0  is  a  much  better  local  variable  than  x,  which  suffers  from  ambiguity 
about  the  location  of  the  "virtual  origin"  of  the  boundary  layer.  Since 
6  was  ignored  in  the  present  analysis,  no  formula  of  this  type  arose. 

It  should  be  pointed  out  again  that  formulas  based  on  R0  cannot  stand 
on  their  own,  because  0(x)  is  not  a  known  geometric  variable  and  must 
be  computed  as  part  of  a  Kaman  integral  type  of  analysis. 

For  computing  the  drag  from  ,  we  have  the  following: 

a)  Blaaius  power-law: 

CD  -  0.074/T^0'2  +  25  % 

b)  Prandt 1 -Schlichting i 

C0  •  (MSS/tlog^)2*5®  «  3  \  l6] 

C)  Schults-Grunow: 

C_  «  0,427/nr»«  r= 
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d)  Karman-Schoenherrj 


I 

3 

i 

J 

.4 

| 

C^11  -  4.13  109^  CD)  ♦  2  %  >  j 

The  Karman/Schoenherr  and  Prandtl/Schlichting  formulas  are  clearly 
superior. 

To  achieve  explicit  formulas  with  ev&n  better  accuracy,  the  writers 
fit  the  numerical  values  from  Eqns.(57)  and  (58)  to  Prandtl/Schlichting 
type  curves,  with  the  following  results? 


a) 

b) 

c) 

d) 


0.225/(log10Rx)2*32  +  0.5  % 

1  64 

0.0253/(logi0RQ)i%OH  +  1*5  % 

0.430/dog^)2*56  +  0.8  % 

O.O385/(logloR0(L)) 1,807  +  0.6  % 


(62) 


These  new  expressions  are  the  most  accurate  simple  and  explicit  formulas 
known  to  the  writers.  They  are  recommended  for  general  usage  for  the 
incompressible  flat  plate  and  for  use  with  the  compressible  flow 

transformations  listed  in  Table  1. _ _ —  - 

The  next  section  will  consider  cases  where  the  freestream  and  wall 


conditions  are  variable. 


V.  COMPRESSIBLE  FLOW  WITH  A  PRESSURE  GRADIENT 


The  chief  use  of  the  present  method  is  in  computation  of  the  skin 
friction  distribution  C^(x)  in  a  turbulent  boundary  layer  with  arbitrary 
distributions  of  M  (x) ,  T  (x) ,  and  T  (x) .  Very  few  existing  methods 

S  6  W 

even  apply  to  such  general  conditions ,  and  these  few  are  all,  to  our 
knowledge,  an  order  of  magnitude  more  complicated  than  the  present 
theory.  The  present  analysis  consists  only  of  a  single  first  order 
differential  equation  for  the  skin  friction,  Eqn.(25),  with  the  various 
coefficients  in  this  equation  being  computed  from  Eqns. (25,28,29,42) . 

Let  us  rewrite  these  equations  here  for  summary  purposes: 


RL  V  -  V’F/V  + 

X4H(1A) '  '/Rl 

dx* 

G  -  3  a  H 

/ 

where: 

a  =  X3(1/V)'Al 

r 

*L 


<u  L/V)  <  u/u)  (T/r  ) 
o  e  e  w  e  w 


5  "  e-  -»<■■&?? 
H  *  °'062 


n  #.  r*  r  •»  n 
- C  — 


(63) 


and  f  4 


1  +  0.11r(k-l)M^(Tft/Tw) 
1  +  0.3(Taw-T„)/Tw 


•f 

Since  the  above  correlations  for  G  and  H  contain  the  thickness  6  , 
we  must  (apologetically)  supplement  Eqn.(25)  with  the  law-of-the-wall 

+  L 

relation,  Eqn.(37),  which  relates  j  to  the  skin  friction  X  ■  (2/Cf)  : 
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X  =  S  +  q  sin  {  0  +  x  t  2(P-Po)  +  ln(|^  |2li)  ]>  , 


where  0  *  sin"^2^  —  -fy  ,  Q  =  ( fiP  +4 y)  ^ ,  and  P  =  (1+  a5+)1*. 


To  match  at  very  low  Reynolds  number  with  the  logarithmic  law -of -the -wall, 

the  initial  conditions  were  taken  to  be  (u+, 6+)  =  (10.0,6.0).  Note  that 

o  o 

1  appears  only  on  the  left  hand  side  and  S'1"  appears  only  in  the  term  P 
on  the  right  hand  side.  However,  in  general,  6  and  Y  are  not  known  in 
advance  and  must  be  computed  from  Eqn.(43)  and  the  local  skin  friction: 


Y  -  r(%i)M2/A2 


(Taw Aw)  -  1 


<TeAw> 


Thus  it  is  definitely  necessary  to  iterate  Eqn. (64)  to  compute  6  ( X) . 

A  FORTRAN-IV  program  is  listed  in  the  Appendix  which  integrates 

Eqn. (63)  subject  to  Eqn. (64,65)  when  the  user  specifies  1)  an  initial 

value  A  at  some  position  x*  ■  x  A;  and  2)  known  distributions  of 
o  o  o 

M  (x) ,  T  (x),  and  T  (x) .  The  program  assumes  a  perfect  gas,  so  that  the 
e  e  w 

velocity  distribution  is  coniputwu  from  U  -  M  (km1  )  ,  and  the  velocity 

6  8  8 

ratio  is  given  by  V  ■  U  /U  ■  (M  /M  )  (T  A  )\  The  computation  of 

e  o  e  eo  e  eo 

Y  from  Eqn, (65)  also  uses  a  perfect  gas  assumption,  and  the  user  would 
be  required  to  modify  these  two  portions  for  real  gas  applications. 

Now  let  us  consider  some  particular  cases. 
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Flow  with  a  Modest  Pressure  Gradient: 

Suppose  that  the  Mach  number  (or  freestream  velocity)  variation  is 
only  slight  and  the  wall  temperature  nearly  constant.  Then  the  parameter 
will  be  very  small  and  we  may  neglect  terms  involving  a  and  o'  in 
Eqn.(63).  We  may  also  forgo  computing  <$+  from  Eqn.  (64) ,  since  it  appears 
only  in  conjunction  with  a  .  Finally,  6  and  y  from  Eqn. (65)  would  be 
roughly  constant.  Equation  (63)  reduces  to: 

12  *  ^-i-5..t53..,G,  V'jfl  ,  (6i 

dx*  G  '  ' 

and  G  =  8.5  exp(0.475  f  X  (T  /T  )**) 

e  w 


We  have  replaced  F  by  (5.53  G).  Since  G  is  approximately  proportional 
to  e*  ,  Eqn. (66)  has  a  closed  form  solution: 


Cf(x*) 


0.42  f2  (Ta/Tw) 

In2 (0.056  f  (T  /T  ) 
e  w 


) 


(MODEST 

PRESSURE 

GRADIENT) 


(67) 


X  A 

where  R  -  (U  L/v  )  V-2,57  f  V+3,57  dx*  . 
err  o  e  1 

0 


Equation  (67)  is  the  compressible  flow  analog  of  the  incompressible 
relation  of  the  same  form  discovered  by  White  (74).  Note  that  it 
merely  modifies  the  incompressible  relation  by  the  same  factors  Fq 
and  fRx  defined  earlier  for  the  flat  plate  in  Eqn. (54) .  Thus  it  is 
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proved  that,  for  modest  pressure  gradients,  the  flat  plate  stretching 
factors  can  be  applied  directly  to  an  incompressible  pressure  gradient 
calculation,  in  the  manner  of  the  Coles  (18)  and  Mager  (44)  compressibility 
transformations.  But  the  idea  breaks  down  entirely  if  the  terms  involving 
a  are  not  negligible.  This  would  explain  why,  as  discussed  by  McDonald 
in  Bertram  (5),  the  simple  integral  transformation  theories  such  as  Pish 
and  McDonald  (26)  and  Reshotko  and  Tucker  (59)  are  accurate  for  modest 
pressure  gradients  but  fail  when  the  gradients  are  strong.  The  authors 
have  found  no  explicit  criterion  for  the  validity  of  Eqn. (67) ,  and  its 
use  in  practical  cases  is  probably  very  limited. 

Flow  with  a  strong  Adverse  Pressure  Gradient: 

If  the  pressure  gradient  is  strong,  we  are  required  to  attack  Eqn.f{63) 
directly  with,  say,  the  computer  program  in  the  Appendix.  To  assess  its 
accuracy,  it  is  desirable  to  have  skin  friction  data  in  a  strong  adverse 
gradient.  The  writers  have  found  only  one  suitable  experiment:  the  flow 
past  a  waisted  body  of  revolution  studied  by  Winter,  Smith,  and  Rotta  (77). 
Although  the  flow  was  axisymmetric,  the  boundary  layer  for  x  greater 
than  24  inches  was  approximately  two-dimensional,  and  we  will  not  consider 
axisymmetric  effects  here.  There  are  other  experiments  -  e.g.  McLafferty 
and  Barber  (46a)  -  which  have  been  compared  with  other  theories.  In 
Bertram  (5),  McDonald  considers  three  such  e>qperiments ,  all  of  which 
measure  only  momentum  thickness  and  shape  factor,  not  skin  friction.  This 
is  almost  unbelievable,  until  we  reflect  again  that  presently  both  theory 
and  experiment  are  locked  in  the  grip  of  the  Karman  integral  relation. 

The  only  possible  reason  a  designer  would  want  to  know  9  or  is  that 
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their  product,  the  displacement  thickness,  is  at  least  nominally  useful, 

to  the  extent  that  it  correlates  such  peripheral  phenomena  as  leading 

edge  shock  wave  interactions  (which  are  not  likely  to  be  turbulent  flow) 

and  local  wall  pressure  fluctuations  (which  correlate  equally  well  with 

the  parameter  6+  computed  in  this  analysis  by  Eqn.(64)).  Nor  is  the 

typical  designer  liable  to  pay  any  more  than  lip  service  to  the  idea  of 

adding  the  displacement  thickness  6* (x)  to  the  body  shape  for  improved 

aerodynamic  computations.  Rather,  the  writers  believe  that  the  only 

parameter  of  primary  design  importance  is  the  skin  friction  Cf(x),  and 

it  seems  incredible  that  an  experiment  could  neglect  this  all  important 
* 

measurement. 

We  consider  now  the  data  of  Winter,  Smith,  and  Rotta  (77).  The 

freestream  distributions  M  (x) ,for  six  different  test  section  Mach 

e 

numbers  (labelled  M  )  are  shown  in  Figure  8.  Since  the  leading  edge 

m 

was  a  thin  cone  and  not  well  approximated  by  the  two-dimensional  equations, 
we  begin  the  computation  at  x  =  24  inches,  where  an  adverse  pressure 
gradient  begins  and  later  levels  out  to  nearly  constant  velocity  at  about 
x  *  45  inches.  The  walls  of  the  model  were  essentially  adiabatic.  The 
curve-fit  velocity  distribution  for  the  present  theory  was  chosen  to  be 
of  the  form 

U.  (x)  *  a  +  (b  +  cx)  e"*3  *  “  ®  X  ,  (68) 

e 

where  (a,b,c,d,e)  were  fitted  constants.  Equation  (63)  was  then  solved 

*Tha  lack  of  friction  data  is  particularly  annoying  in  view  of  the  fact 
that  the  "measurement"  of  8  and  H]_2  actually  involves  a  complete 
survey  and  integration  of  both  the  velocity  and  temperature  profiles 
across  the  entire  boundary  layer  at  each  station. 
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X  -  in. 


Figure  8.  FREESTREAM  MACH  NUMBER  DISTRIBUTIONS 
ON  A  WAISTED  BODY  OF  REVOLUTION,  FROM 
THE  DATA  OF  WINTER,  SMITH,  AND  ROTTA  (77). 


for  Cf (x)  on  the  IBM  360/50  digital  computer  at  the  University  of  Rhode 
Island.  The  initial  condition  was  taken  to  be  the  measured  skin  friction 
at  x  =  24  inches.  A  complete  run  for  a  given  on  the  computer  took 
about  ten  seconds,  the  limitation  being  time  required  for  print -out. 

The  comparison  between  theory  and  experiment  is  shown  in  Figures  9  and 
10.  Also  shown  are  the  finite  difference  confutations  of  Herring  and 
Mellor  (30).  The  present  theory  is  seen  to  be  in  good  agreement  and  in 
fact  is  superior  in  every  case  to  the  much  more  complex  analysis  of 
Herring  and  Mellor  (30).  Like  most  other  methods,  Herring  and  Mellor 
key  their  initial  condition  to  the  measured  momentum  thickness  and  shape 
factor  and  use  these  two  to  confute  the  initial  skin  friction,  which 
happened  to  fall  much  too  high  at  the  larger  Mach  numbers.  The  same 
difficulty  was  reported  in  the  integral  method  of  Alber  and  Goats  (2)  and 
in  the  review  paper  by  McDonald  in  Bertram  (5).  The  present  theory,  of 
course,  keys  directly  to  the  skin  friction  -  a  decided  advantage  over 
Karman-oriented  methods.  In  no  case  did  the  present  theory  predict 
separation,  although  the  highest  Mach  number  run  (M^1  2.8)  hinted  of  a 
near-separation  condition  with  an  initial  decrease  in  the  denominator 
(G  -  3  a  H)  of  Eqn. (63) .  The  agreement  of  the  present  theory  in  the 
relaxation  sons  at  the  trailing  edge  is  surprisingly  good,  considering 
that  relaxation  is  the  chief  area  of  inaccuracy  in  the  related  incompressible 
analysis  of  white  (74).  Perhaps  relaxation  is  not  as  serious  a  problem 
for  a  law-of -the -wall  approximation  in  supersonic  boundary  layers. 
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Figure  9.  COMPARISON  OF  THEORY  WITH  THE  DATA 
OF  WINTER,  SMITH,  AND  ROTTA  (77).  Part  1. 


Flow  with  a  Strong  Favorable  Pressure  Gradients 

Skin  friction  data  in  a  supersonic  favorable  gradient  were  recently 

reported  by  Brott  et  al  (9).  Using  a  flexible  nozzle,  these  workers 

generated  a  freestream  which  increased  from  M  =  3  to  about  M  *  4.6  1 

e  e  3 

in  a  distance  of  sixty  inches.  The  measured  freestream  Mach  numbers  are  : 

1  /*»  •  f 

shown  in  Figure  11-a.  They  roughly  approximate  a  power-law  Mg  =  1.0?  x  , 

which  was  used  as  a  curve -fit  in  the  present  theory,  Eqn.  (63),  to  compute 

Cf(x),  beginning  at  x  =  44  inches.  The  data  was  taken  for  a  range  of  values  ' 

of  the  tunnel  stagnation  pressure,  which  effectively  corresponds  to  a  family 

of  values  of  the  nominal  Reynolds  number  because  of  the  variation  in 

freestream  density.  Figure  11-b  compares  the:  theoretical  and  experimental 

(wall  shear  balance)  measurements  of  skin  friction  for  four  stagnation 

pressures.  Also  shown  is  a  flab  plate  computation  for  pQ  =  150  psia, 

which  illustrates  the  usual  fact  that  favorable  gradient  friction  lies 

above  the  equivalent  flat  plate  values.  The  wall  temperature  was  slightly 

cold,  T  =  0.82  T  ,  and  the  theory  was  run  for  this  condition.  The 
w  aw  ■* 

i 

agreement  of  the  theory  is  good  except  that  it  falls  somewhat  low  at  the 
trailing  edge.  Friction  data  by  Lee  et  al  (42)  at  zero  pressure  gradient 
in  the  same  wind  tunnel  also  rise  somewhat  higher  downstream  than  a  flat 
plate  calculation.  The  reason  for  this  slight  discrepancy  is  not  known 
to  the  writers. 

An  interesting  approximation  for  strong  favorable  gradients  at  large 

•  "  '  ..  •  j- 

Reynolds  was  found  from  an  inspection  of  the  computer  results.  If  the 
gradient  is  truly  strong  (large  negative  a),  the  terms  involving  H  are 

j 

•  •  .  -  j 

dominant  in  Eqn.  (63),  which  simplifies  to:  I 

1 

j 
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Mach  number. 


b)  Local  Skin  Friction. 


Figure  11.  COMPARISON  OF  THE  PRESENT  THEORY  WITH  THE 
FAVORABLE  GRADIENT  DATA  OF  BROTT  ET  AL  (9). 


dX 

dx* 


(69) 


x4  (1/V)»*/Rl  _  -  X  (1/V) l » 

-  3  o  ~  3  (1/V)  * 


This  equation  integrates  exactly  into  the  following  simple  approximation: 


(STRONG 

FAVORABLE 

GRADIENT) 


(70) 


In  this  simple  theory,  the  skin  friction  depends  only  upon  the  distribution 

of  (U  /U  ) '  and  is  unaffected  by  the  level  of  Mach  number,  Reynolds 
o  e 

number,  or  wall  temperature.  Equation  (70)  is  also  shown  in  rigure  11-b, 
and  the  agreement  is  seen  to  be  best  at  high  Reynolds  nuthers. 

Some  Idealised  Sample  Calculations: 

The  experiments  discussed  above  did  net  indicate  flow  separation  or 
even  suggest  trends  for  theoretical  comf  ;»rison.  It  was  decided  to  complete 
this  section  by  developing  a  fer.:  idealized  cases  for  which  the  solution 
of  the  present  theory,  Eqn. (63),  could  predict  the  effect  of  certain  flow 
parameters  on  adverse  pressure  gradients.  The  rase  selected  was  an 
adiabatic  wall  with  an  exponential  frea&tream  Mach  number  distribution: 

M  (x)  ='  "  M  ,  (71) 

e  o 

and  computer  runs  were  made  for  various  values  of  Mq,  k,  and  -he  Reynolds 
number  (which  can  be  interpreted  from  Eqn. (26)  as  either  actual  Reynolds 
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number  change  or  as  a  change  in  the  effective  wall  temperature  ratio.  The 
results  of  these  idealised  computations  are  shown  in  Figures  12,  13,  and 
14,  showing  the  effect  of  K,  Mq,  and  R^  respectively. 

Figure  12  shows  that,  for  a  given  Mach  number  and  Reynolds  number, 
an  increase  in  the  adverse  gradient  (K)  will  eventually  cause  flow 
separation  (Cf  =  0).  Confutations  at.  other  values  of  Mq  and  R^  confirm 
this  effect.  It  appears  that  a  suitably  relentless  adverse  gradient  will 
ultimately  drive  any  turbulent  boundary  layer  to  separation. 

Figure  13  shows  the  effect  of  Mach  number  level  on  these  idealized 
flows.  An  increase  in  Mach  number  naturally  tends  to  drive  the  skin 
friction  to  lower  values,  just  as  in  the  flat  plate  case,  but  no  tendency 
toward  separation  is  noted  until  Mq  reaches  hypersonic  values  of  the 
order  of  ten.  This  effect  is  partly  computational  in  nature,  because  the 
present  theory  predicts  the  onset  of  separation  at  a  finite  value  of  the 
order  of  0.0001,  which  is  more  likely  to  happen  for  =  10.  we  may 
speculate,  however,  that,  for  given  gradient,  an  increase  in  Mach  number 
appears  to  increase  the  tendency  toward  flow  separation. 

Finally,  Figure  14  shows  the  effect  of  Reynolds  number  for  a  given 
Mach  number  and  gradient.  No  flow  separation  is  indicated  at  any  level 
of  R^.  if  anything,  the  "tendency"  toward  flow  separation  is  stronger 
at  the  low  Reynolds  numbers.  This  was  also  the  case  in  the  incompressible 
flow  analysis  of  White  (74) . 

It  is  hoped  that  more  skin  friction  data  under  variable  Mach  number 
and  wall  temperature  conditions  will  be  forthcoming  from  workers  in  the 
field  of  compressible  turbulent  boundary  layers. 
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Figure  12.  THEORETICAL  EFFECT  OF  PRESSURE  GRADIENT 
MAGNITUDE  ON  AN  IDEALIZED  RETARDED 
SUPERSONIC  BOUNDARY  LAYER,  FROM  EQN.  ( 63) . 
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Figure  13.  THEORETICAL  EFFECT  OF  MACH  NUMBER 
MAGNITUDE  ON  AN  IDEALIZED  RETARDED 
SUPERSONIC  BOUNDARY  LAYER,  FROM  Eqn.  (63) 


x/L 

Figure  14.  THEORETICAL  EFFECT  OF  REYNOLDS  NUMBER 
MAGNITUDE  ON  AN  IDEALIZED  RETARDED 

SUPERSONIC  BOUNDARY  LAYER,  FROM  EQN.(63,. 
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VI.  CONCLUSIONS 


It  is  the  purpose  of  this  report  to  develop  an  entirely  new 
approach  to  the  calculation  of  skin  friction  in  compressible,  two- 
dimensional  turbulent  boundary  layers  with  arbitrary  pressure  gradients 
and  arbitrary  heat  transfer.  The  analysis  uses  the  familiar  law-of-the- 
wall,  u^(y+),  which  is  generalized  in  Section  III  to  include  the  effect 
of  local  pressure  gradient,  heat  transfer,  and  compressibility.  The 
final  expression  for  the  chosen  form  of  the  law-of-the^wall  is  given  by 
Eqn.  (37) .  By  utilizing  this  expression  as  an  "equation  of  state"  of 
turbulence,  the  momentum  and  continuity  equations  can  be  integrated 
across  the  boundary  in  law-of -the -wall  coordinates.  The  result  is  a 
single  differential  equation  with  the  skin  friction  Cf(x)  as  the  only 
variable,  Eqn.  (63).  This  is  in  sharp  contrast  to  the  well  known  Karman 
integral  approach,  which,  by  integrating  in  physical  rather  than  wall- 
related  coordinates,  introduces  not  only  skin  friction  but  also  other 
(unwanted)  variables  such  as  integral  thicknesses  and  shape  factors. 

The  present  theory  is  much  simpler  than  any  Karman-type  theory  and  is 
frankly  meant  to  compete  with  or  eventually  replace  the  Karman  approach. 
A  computer  program  for  the  new  theory  is  given  in  the  Appendix,  but 
the  computations  are  simple  enough  to  be  performed  by  hand. 

Some  of  the  highlights  of  the  new  theory  are  as  follows t 

1.  It  is  valid  (but  approximate)  for  arbitrary  freestream  and 
wall  temperature  distributions. 

2.  The  skin  friction  is  the  only  variable,  and  the  momentum 
thickness,  the  displacement  thickness,  the  shape  factor,  etc.  are  not 
required .  Nor  is  any  skin  friction  correlation  or  cosqpressibility 
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transformation  required. 

3.  The  analysis  contains  an  explicit  separation  criterion. 

4.  For  flat  plate  flow,  the  analysis  leads  to  a  simple  algebraic 
expression  for  skin  friction  which,  for  adiabatic  walls,  has  the  smallest 
mean  absolute  error  of  any  know  flat  plate  theory.  For  a  flat  plate 
with  cold  walls,  the  present  theory  takes  a  creditable  fourth  place 
behind  the  theories  of  Moore  (52),  Spalding  and  Chi  (66),  and  Van 
Driest  (72).  For  heat  transfer,  then,  the  chief  drawback  of  the  present 
theory  is  the  lack  of  an  adequate  temperature  law-of -the -wall. 

5.  For  compressible  flow  with  very  modest  pressure  gradients, 
the  present  theory  leads  to  a  closed  form  solution  for  C^(x)  which 
gives  insight  into  why  some  other  well  known  theories  are  accurate  for 
mild  pressure  gradients  but  fail  when  gradients  are  strong. 

6.  For  compressible  flow  with  strong  pressure  gradients,  '"-he  new 
theory,  Eqn.  (63) ,  apparently  gives  the  best  agreement  with  experiment 
of  any  known  compressible  turbulent  boundary  layer  theory.  However, 
there  are  so  few  data  that  to  term  the  new  theory  "most  accurate"  is 
obviously  a  premature  move. 

Overall,  it  is  felt  that  the  new  theory  is  definitely  useful  from 
the  point  of  view  of  1)  simplicity;  2)  accuracy;  and  3)  insight  into 
the  behavior  of  turbulent  skin  friction  in  high  speed  flows.  It  is 
also  felt  that  this  new  approach  can  be  readily  extended  to  include 
more  accurate  heat  transfer  computations  and  other  geometries. 
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APPENDIX  I . 

A  FORTRAN  PROGRAM  FOR  SOLVING  EQUATION  (63) 


C  FRANK  M.  WHITE  -  GEORGE  H.  CHRISTOPH  -  UNIVERSITY  OF  RHODE  ISIA1ID 
C  THIS  PROGRAM  SOLVES  EQUATION  (63)  WITH  A  RUNGE-KTJTTA  SUBROUTINE. 

REAL  M,  N,  MZ 
DIMENSION  Y (30) ,  Z(30) 

ASIN(X)  =  ATAN(X/SQRT(1.-X‘X)) 

C  READ  IN  VALUES  OF  UZERO,  MZERO,  TZERO,  AND  CURVE-FIT  CONSTANTS. 

6  READ ( 5,1)  UZ,  MZ,  TZ,  A,  B,  C,  D 

1  FORMAT (7F10. 2) 

L  =  0 

C  READ  IN  RL,  XZERO,  XMAX,  DELTA-X,  MU  EXPONENT,  DELTA PLUS,  LAMBDA-ZERO. 
READ (5,1)  R,  X,  XMAX,  DX,  N,  DEL,  Y(l) 

WRITE(6,61)  MZ,  Y(l) 

61  FORMAT ( ' 1  THIS  RUN  IS  FOR  MZERO  -  ' ,F9.3,9X, 'AND  LAMBDA  *\F9.2/) 

C  TEST  AND  ENTER  SUBROUTINE  IF  XMAX  HAS  NOT  BEEN  EXCEEDED. 

47  IF(X-XMAX)  2,  16,  16 

2  CALL  RUNGE (1 ,Y,Z,X,DX,L,K) 

GO  TO  (10,20) ,  K 

C  NOW  COMPUTE  THE  VELOCITY  V  AND  ITS  TWO  DERIVATIVES  FROM  CURVE -FITS. 

10  V  *  A  +  B*X  +  C*X*X 
VP  »  B  +  2 .*C*X 
VPP  -  2.‘C 

C  COMPUTE  THE  TEMPERATURE  OF  WALL  AND  FRBESTRBAM  FROM  CURVE-FITS. 

TE  »  B  +  C*X  +  D*X“3 
TW  -  A  +  B*X‘X 

C  COMPUTE  THE  FIRST  AND  SECOND  DERIVATIVES  OF  (1/V) . 

UP  *  -VP/V/V 

UPP  -  (-VPP  +  2.‘VP*VP/V)/V/V 
C  COMPUTE  MACH  NUMBER,  ALPHA,  BETA,  GAMA,  AND  UPLUS-E. 

M  -  MZ*V*SQRT(TZ/TE) 

RL  «  R*(TE/TW)**(.5-HI) 

ALF  -  Y(1)“3*UP/RL  +  .0000001 

BETA  -  ((TE/rW)*(l.+.l78*M*M)  -  1.) /SQRT(TE/TW) 

GAMMA  "  .  178‘M*M/Y(1)“2 
UPLUS  -  Y(l)*SQRT(TE/rW) 

C  NOW  USE  THE  LAW-OF -THE -WALL,  EQN(37) ,  TO  COMPUTE  DBLTAPLUS  BY  ITERATION. 
PZ  ■  SQRT(1.  +  6. *ALF) 

PR  »  (PZ  -l.)/(PZ  +  1.) 

Q  -  SQRT( BETA* BETA  +  4. ‘GAMMA) 

SZ  -  ASIN( (20. ‘GAMMA  -  BETA) A}) 

IF(UPLUS  -  (Q+B8TA)/2. /GAMMA)  41,  42,  42 

42  S  «  1.570796 
GO  TO  43 

41  S  -  AS1N( (2 «‘GAMMA‘UPLUS  -  BETA)  A}) 

43  8  •  .4* (S-SZ) /SQRT(GAMMA) 

IP(A3S (ALF) -.00001)  40,  40,  81 

40  DEL  -  6.*EXP(S) 

GO  TO  4 

81  P  *  SQRT(ABS(1.  +  DEL* ALF)) 
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FORTRAN  PROGRAM  FOR  EQUATION  (63)  -  Continued 


DO  80  I  =  1,5 
T  =  PR»EXP(S-2.*(P-PZ)) 

P  *  P  -  ((P-l)/(P+l)-T)/(2./(P+l)**2+2.*T) 

IF(P)  26,  26,  80 
80  CONTINUE 

DEL  *  (P*P-1)/ALF 
IF(ALF*DEL  +  1)  26,  4,  4 
26  DEL  «  -1/ALF 
4  AD  *  ALF*DEL  +  .0000000001 

C  NOW  COMPUTE  THE  FUNCTIONS  G,  F,  AND  H  FOR  EQN. (63). 

FAC  =  (l.+.22*GAMMA*UPLUS**2)/(l.+.3*BETA*UPLUS) 

G  »  8 . 5*EXP ( . 475*FAC*UPLUS/ ( 1 .  +.  1*AD/ABS (AD) ** . 5) ) 

H  *  .062*EXP(.84*FAC*UPLUS/(1.+.12*AD/ABS(AD)**.4)) 

F  *  5.53*G 

C  FINALLY,  COMPUTE  THE  DERIVATIVE  OF  LAMBDA  FROM  EQN. (63) . 

GNET  -  G  -  3.*ALF*H 

2(1)  -  (RL*V  -  VP*FA  +  Y(1)**4*H*UPP/RL)/GNBT 
C  RETURN  TO  THE  SUBROUTINE  UNLESS  GNET  HAS  BECOME  NEGATIVE  (SEPARATION) . 
IF (GNET)  606,  606,  2 

C  THE  SECOND  SUBROUTINE  BRANCH  POINT  (20)  PRINTS  OUT  THE  ANSWERS. 

20  CF  -  2./Y(l)*#2 

WRITE (6, 361  X,  DEL,  CF,  V,  GNET 

36  FORMATC  X*' ,F9.3,  ’  DEL«\P9.1,'  CF•,,F9.6,  ,  V-*,F9.5,'  GN-\F9.0) 
GO  TO  47 
16  WRITE (6,86) 

86  FORMAT ( '0  NO  SEPARATION  DURING  THIS  RUN. ') 

GO  TO  6 

606  WRITE (6.96)  X 

96  FORMAT ( '0  SEPARATION  HAS  OCCURRED  AT  X  -  ',F10.5) 

GO  TO  6 
END 


NOTEt  This  program  should  be  coafeinsd  with  the  Subroutine  ROUGE 
which  is  printed  end  discussed  on  the  next  two  pages. 
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RUNGE-KUTTA  SUBROUTINE 

For  convenience  to  the  reader,  a  description  of  the  arguments 
of  the  Runge  Kutta  subroutine  and  the  subroutine  itself  are  given 
below. 

CALL  RUNGE  (  N,  Y,  F,  X,  H,  M,  X) 

N  -  is  the  number  of  differential  equations  to  be  solved  (set  by 
the  programmer) 

Y  -  is  an  array  of  "N"  dependent  variables  (initial  values  set  by 
the  prograwaer) 

F  -  is  an  array  of  sise  "N"  of  the  derivatives  of  the  variables  "Y" 
(expression  for  each  F(l)  must  be  given  by  the  programmer) 

X  -  is  the  independent  variable  (initialised  by  the  programmer) 

H  -  is  the  desired  increment  in  X  (given  by  the  programmer) 

M  -  is  an  index  used  in  the  SUBROUTINE  which  must  be  set  equal  to 
sero  by  the  progreemer  before  the  first  CALL 
K  -  is  an  integer  calculated  by  the  SUBROUTINE  which  is  used  in  a 
computed  GO  TO  statasmnt  in  the  main  program.  For  example,  the 
pxogxaMMr  may  use  GO  TO  (10,  20) ,  K  where  statement  10  calculates 
the  derivatives  F(I)  and  statement  20  prints  the  answers  X  and  Y(I) . 
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SUBROUTINE  RUNGE  (N,Y,F,X,M,K) 

THIS  RUNGE-KUTTA  SUBROUTINE  WAS  DEVELOPED  BY  F.M.  WHITE 
AND  PERFORMS  RUNGE-KUTTA  CALCULATIONS  BY  GILLS  METHOD 
DIMENSION  Y(10) ,  F(10),  Q(10) 

M  -  M  +  1 

CO  TC  (1,4, 5, 3, 7) ,  M 

1  DO  2  T  -  1,  N 

2  Q(I)  =  0. 

A  =  .5 
GO  TO  9 

3  A  -  1.707107 

IF  YOU  NEED  MORE  ACCURACY,  USE  A  =  1.7071067811865475244 

4  x  =  x  +  .5*h 

5  DO  6  I  =  1,  N 

Y(I)  •  Y(I)  +  A*(F(I)*H-Q(I)) 

6  Q(I)  =  2.*A*H*F(I)  t  (1.  —  3.*A)*Q(I) 

A  «  .2928932 

IF  YOU  NEED  MORE  ACCURACY,  SET  A  =  .2923932136134524756 
GO  TO  9 

7  DO  8  I  =1,  F 

8  Y(I)  *  Y(I)  +  H*F(lj/6.  -  Q(I)/3. 

M  »  0 

K  =  2 
GO  TO  10 

9  K  =  1 
10  RETURN 

END 
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APPENDIX 


III 


EXPERIMENTAL  DATA  CP  SKIN  FRICTION  COEFFICIENTSAT 
REYNOLDS  NUMBERS  AND  MACH  NUMBERS 
(For  Adiabatic  Flow) 


M 

e 

*0 

R  X  10’6 

X 

Cf  (exp) 

Brinich  and  Dlaconis 

(8) 

3.05 

2130 

.00206 

3.05 

2600 

.00199 

3.05 

3030 

.00190 

3.05 

4400 

.00174 

3.05 

7800 

.00146 

3.05 

10300 

.00142 

3.05 

13000 

.00129 

3.05 

15400 

.00128 

Chapman  and  Kester 

(15) 

.91 

4704 

.81 

4.82 

.81 

6.20 

.81 

6.68 

.81 

7.42 

.81 

7.78 

.81 

8.20 

.81 

9.00 

.81 

9.80 

.81 

10.90 

.81 

11.90 

.81 

T'.r  .  - 

.  12.00 

-  ;,j; --  ■ 

e81 

13,25 

.61 

13.50 

.81 

■  : 

15.00 

.81 

15.4 

.81 

17.6 

.81 

17.3 

.81 

....  -■  • 

18.0  '  >' 

■ 

.81 

20.7 

.61 

..  .  -  ' 

23,3 

.81 

26.9 

.81 

31.8 

2.5 

5.78 

2.5 

6.98 

2.5 

7.7 

2.5 

9.0 

2.5 

9.3 

2.5 

11.2 

2.5 

-  _ 

11.3 

VARIOUS 


Cj  (exp) 


.00308 

.00293 

.00280 

.00268 

.00282 

.00279 

.00259 

.00276 

.00268 

.00261 

.00263 

.00264 

.00260 

.00261 

.00256 

.00250 

.00250 

.00251 

.00250 

.00243 

.00240 

.00237 

.00232 

.00216 

.00202 

.00202 

.00195 

.00195 

.00191 

.00189 


2.5 

2.5 

2.5 

2.5 

2.5 

2.5 

2.5 

2.5 

2.5 

2.5 

2.5 

2.5 

3.6 
3.6 
3.6 
3.6 
3.6 
3.6 


1.966 

1.978 

1.982 

2.540 

2.568 

2.578 

3.690 

3.701 

3.697 

4,512 

4.554 

4.545 

4.504 

4.544 


2.5 

2.5 

2.5 


.352 

.630 

.764 

1.24 

1.26 

1.37 

1.44 


>  Rx  X  10~6  Cf  <exP) 

Chapman  and  Kcster  (cont . ) 

12.4 

14.3 
1414 
16.0 
16.0 
17.9 
18.0 
21.0 

24.2 

26.4 

28.3 

31.2 
6.25 
9.46 
10.0 

16.2 

17.5 
18.3 

Coles  (18) 

.00272 
.00218 
.00202 
.00242 
.00181 
.00166 
.00211 
.00162 
.00138 
.00148 
.00122 
.00131 
.00155 
.00126 

Cope  (20) 

.00250 
.00210 
.00210 

Dhawan  (22) 


2980 

6470 

8570 

2190 

6600 

10200 

2120 

4100 

7560 

3470 

6590 

4980 

2900 

5240 


2565 

2720 

5655 


Cd  (exp) 


.00190 

.00190 

.00188 

.00182 

.00182 

.00183 

.00179 

.00174 

.00170 

.00168 

.00163 

.00167 

.00170 

.00163 

.00165 

.00153 

.00154 

.00151 


1 

1 

1 

1 

1 

1 

i 


.00346 

.00340 

.00329 

.00296 

.00296 

.00282 

.00279 
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M 

X. 


Rx  X  10~6  Cf  <CXp)  Cd  (exp) 

phawan  (cont.) 

1  .00272 


0.7 

0.7 

0.7 

0.7 

3.07 

3.07 

3.07 

3.07 

3.07 

3.7 

3.7 

3.7 

4.54 

4.54 

4.54 

4.54 

4.54 

4.54 

4.54 

4.54 

4.54 

4.54 

4.54 

4.54 

4.54 

4.54 

4.54 


.18 

.20 

.26 

.27 

.31 

.33 

.37 

.46 

.49 

.56 

.37 

.65 

.75 

.85 

.97 

1.45 

1.48 


Goddard  (27) 

1.89 
3.80 
5.71 
7.21 

2.54 
2.79 

3.41 

4.30 

5.30 
3.68 

4.40 
5,05 
4.23 

4.42 
4.45 

4.55 

4.90 

5.20 

5.40 

6.20 

6.40 

6.40 
7.10 

7.40 
7.40 

7.90 

8.20 

HalcHinen  (28) 


.33 

.00417 

.36 

.00409 

.48 

.00395 

.48 

.00380 

.55 

.00384 

.59 

.00374 

.63 

.00366 

.76 

.00340 

.82 

.00344 

.90 

.00336 

.90 

.00337 

1.00 

.00330 

2.03 

.00301 

1.12 

,.00301 

1.20 

•00300 

1.04 

.00300 

1.04 

.00294 

.00412 

.00340 

.00327 

.00328 

.00193 

.00189 

.00182 

.00172 

.00169 

.00169 

.00158 

.00150 

.00150 

.00140 

.00135 

.00149 

.00142 

.00140 

.00139 

.00140 

.00133 

.00129 

.00140 

.00130 

.00123 

.00127 

.00130 
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M 

e 


1.49 

1.50 
1.50 
1.50 
1.50 
1.50 
1.52 
1.52 
1.52 
1.52 

1.71 

1.72 

1.73 

1.73 

1.74 

1.74 

1.75 

1.76 
1.76 


2.445 

2.961 

3.443 

2.468 

2.978 


1.604 
1.592 
2.182 
2.179 
2.188 
2.187 
2.182 
2.185 
*  2.146 
1.595 
1.595 

1.590 

1.591 
1.588 
1.593 
2.115 
2.172 
2.178 
2.192 
2.198 


Rx  x  10  6 
Hakkinen  (cont. ) 


Cf  (exp) 


1.04 

.00300 

1.02 

.00291 

1.04 

.00302 

1.02 

.00302 

1.03 

.00302 

1.03 

.00292 

1.02 

.00309 

1.02 

.00307 

1.01 

.00309 

1.01 

.00307 

.68 

.00324 

.67 

.00321 

.66 

.00321 

.84 

.00313 

.67 

.00323 

.85 

.00319 

.67 

.00323 

.84 

.00310 

.84 

.00312 

Hopkins  and  Keener  (33) 


59680 

55590 

53740 

75260 

68140 


.00126 

.00111 

.00091 

.00127 

.C'UO 


Jackson,  Caarmecki  and  Manta  (36> 


80156 
10845 
50989 
43716 
36860 
29198 
20974 
7  132 
7556 
•3872 
72030 
45061 
30511 
17090 
58977 
9657 
13799 
25216 
35099 
44303 


.00162 

.00217 

.001444 

.001461 

.001495 

.001530 

.001614 

•001766 

.001865 

.001620 

.001660 

.001760 

.001860 

.002080 

.001700 

.001654 

.001636 

.001505 

001449 

001400 


Cd  (exp) 
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Jackson  et.  al.  (con».) 


2.200 

52405 

.001404 

2.202 

60112 

.001387 

2.172 

14510 

.001636 

2.159 

35508 

.001449 

2.188 

44672 

.001400 

2.192 

52248 

.001404 

2.192 

59044 

.001387 

2.163 

26336 

.001505 

2.161 

13586 

.001636 

2.110 

9723 

.001654 

2,186 

37982 

.001365 

2.194 

94481 

.001256 

2.188 

63929 

.001289 

2.195 

85542 

.001256 

2.192 

73696 

.001286 

2.189 

61403 

.001289 

2.182 

49008 

.001324 

2.142 

21061 

.001484 

2.083 

13360 

.001544 

2.172 

60685 

.001369 

2.165 

40131 

.001454 

2.138 

17131 

.001690 

2.170 

68974 

.001271 

2.169 

45284 

.001324 

2.115 

19159 

.001522 

2.199 

65500 

.001289 

2.149 

21008 

.001484 

1.587 

123836 

.001408 

1.594 

107318 

1.594 

91405 

.001463 

1.587 

72236 

.001520 

1.575 

50908 

.001623 

1.548 

26317 

.001821 

1.469 

15675 

.001993 

1.599 

78918 

.001550 

1.602 

93874 

.001530 

1.598 

67789 

.001580 

1.593 

52482 

.001650 

1.586 

35864 

.001770 

1.567 

18258 

.001950 

1.555 

11703 

.002080 

1.598 

•9760 

.001530 

1.597 

60359 

*.001616 

1.579 

20816 

.001934 

1.596 

101828 

*001540 

1.596 

67108 

.001590 

1.579 

23077 

Korkagl  (40) 

.001860 

5.787 

2477 

.001316 

75 


001316 


M 


e 


5.770 

5.792 

5.805 


2.95 

2.95 

2.95 

2.95 

2.95 

•2.93 

2.95 

2.95 

2.95 

2.95 

2.95 

2.95 

2.95 

2.95 

2.95 

2.95 

2.95 

2.95 

2.95 

2.95 

4.2 

4.2 

4.2 

4.2 

4.-2 

4.2 

4.2 

4.2 

4.2 

4.2 

4.2 

4.2 

4.2 

4.2 

4.2 

4.2 

4.2 

4.2 

4.2 

4.2 

4.2 

4.2 

4.2 

4.2 

4.2 


\ 


Cf  (exp) 


Korkcgl  (c°nt. ) 


2780 

3429 

4040 


.001275 

.001223 

.001179 


jk?A*An9*  Chapman,  Nyholm  and  Thomas  (46) 


6.18 

.00160 

8.30 

.00155 

9.01 

.00154 

10.50 

.00150 

12.00 

.00146 

13.00 

.00145 

14.8 

.00142 

17.2 

.00140 

19.9 

.00136 

21.9 

.00134 

24.6 

.00133 

26.0 

.00130 

27.5 

.00130 

30.5 

.00130 

34.0 

.00126 

35.5 

.00129 

37.5 

.00124 

42.0 

.00123 

54.0 

.00120 

65.0 

.00190 

4.63 

.00132 

5.90 

.00130 

6.64 

.00126 

7.53 

.0012S 

6.0 

.00120 

9.12 

.00118 

9.8 

.00115 

11.2 

.00114 

12.6 

.00110 

14.0 

.00109 

15.4 

.00107 

17.5 

•00105 

19.5 

.00104 

22.5 

.00101 

25.0 

.000995 

27.5 

.000975 

20.2 

.000960 

24.8 

.000948 

27.6 

.000928 

42.0 

.000922 

48.0 

.000920 

55.0 

.000898. 

61.0 

.000380 

68.2 

.000860 

75.0 

.000855 

Cd  (exp) 
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f 

1 

i 

$■ 

! 

M 

e 

*0 

R  x  10*6 

X 

Cf  (exp) 

Cd  (exp) 

! 

f 

Matting  et.  al.  (cont.) 

l 

4.2 

84.0 

.000840 

4.2 

95.2 

.000828 

6.7 

7.17 

.000680 

6.7 

7.64 

.000660 

6.7 

9.60 

.000634 

T' 

6.7 

14.80 

.000610 

6.7 

18.60 

.000600 

> 

6.7 

27.8 

.000530 

V 

6.7 

33.0 

.000512 

6.7 

45.0 

.000510 

t 

6.7 

46.8 

.000529 

9.9 

14.7 

.000331 

9.9 

20.4 

.000325 

Monaghan  and  Johnson 

(51) 

2.43 

1072 

.596 

.00295 

2.43 

1307 

.772 

.00280 

2.43 

1451 

.824 

.00272 

! 

2.43 

1792 

1.187 

.00250 

2.43 

2521 

1.93 

.00237 

2.43 

3268 

2.36 

.00227 

2.43 

3950 

2.93 

.00218 

Monaghan  and  Cooka 

(50) 

2.82 

886 

.490 

.00272 

2.82 

1493 

.945 

.00243 

2.82 

1777 

1.346 

.00219 

2.82 

2180 

1.570 

.00203 

2.82 

2660 

1.860 

.00196 

' 

.  ; , ' 

2.82 

2789 

2.178 

.00193 

:■■  ■  ■ 

Moor*  and  Karknasi 

(53) 

2.831 

•27 

.000900 

2.843 

611 

.000953 

2.865 

376 

.000987 

1  •• 

2.897 

318 

.001020 

i 

f 

2.787 

1120 

.000849 

' 

} 

2.809 

•73 

.000874 

j 

2.828 

512 

.00094o 

1 

2.669 

1410 

•000862 

; 

2.693 

1180 

.000891 

f 

O'Donall  (56) 

aMMM»*«MMk 

i 

2.41 

1530 

.00240 

j.! 

2.41 

2140 

.00240 

! 

$! 

2.41 

2200 

.00223 

■ii 

2.41 

3000 

.00236 

j> 
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b 

1 

l: 

l 

M 


e 


2.41 

2.41 

2.41 

2.41 

2.41 

2.41 

2.41 

2.41 

2.41 

2.41 

2.41 


2.5 

2.5 

2.5 

2.5 

2.5 

2.5 

2.5 

2.5 

2.5 

2.5 

2.5 

2.5 

2.5 

2.5 

2.5 

2.5 

2.5 

2.5 

2.5 


1.724 

1.724 

1.782 

1.726 

2.017 

1*996 

2.000 

2.249 

2.242 
2.236 

2.243 
2.502 
2.533 
2.541 


R©  Rx  *  N>~6  Cf  (exp) 
0 ’Done 11  (cont.) 

3600  .00190 
3820  .00200 
4100  .00180 


.00173 

.59 

.93 

1.42 

2.12 

2.41 

2.90 

3.26 


Rubcsin,  Maydew  and  Varga  (60) 


6082 

6093 

11644 

19833 

6113 

11015 

20090 

6182 

8301 

10711 

20490 

6085 

9639 

18811 


3.14 

3.49 

3.63 

4.13 

4.13 

4.73 

4.73 

5.45 

6.09 

2.70 

3.21 
3.56 
3.69 
4.20 

4.22 
4.80 
4.83 
5.48 
6.17 

tHutti  (64) 


•002225 

.002225 

.001947 

.001784 

.002060 

.001810 

.001635 

.001985 

.001816 

.001704 

.001623 

.001804 

•001583 

.001560 


C  (exp) 
a 


.00360 

.00326 

.00300 

.00292 

.00291 

.00268 

.00260 


.00266 

.00255 

.00255 

.00246 

.00253 

.00241 

.00242 

.00241 

.00231 

.00247 

.00264 

.00250 

.00251 

.00244 

.00250 

.00241 

.00238 

.00239 

.00229 
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Cd  (exp) 


M 

R 

c 

e 

2.8 

9600 

2.8 

9300 

2.8 

10800 

2.8 

11700 

2.8 

12000 

2.8 

12300 

2.8 

13400 

2.8 

13600 

2.8 

2.8 

2.8 

2.8 

2.8 

2.8 

2.8 

2.8 

1.739 

12490 

1.744 

12240 

1.744 

8429 

1.739 

3589 

1.737 

3443 

2.019 

12320 

2.009 

7528 

2.007 

2899 

2.238 

11670 

2.227 

6892 

2.230 

2520 

2.490 

11400 

2.483 

6097 

2.502 

6072 

2.484 

2660 

2.739 

11900 

2.724 

6304 

2.729 

3048 

2.949 

11400 

2.949 

6041 

2.958 

2740 

3.161 

11310 

3.168 

7149 

3.389 

11270 

3.402 

7785 

3.400 

2758 

3.681 

10180 

3.681 

9836 

3.667 

7991 

3.672 

2096 

3.684 

2075 

*«  *  10'6 


C£  (exp) 


Spivak  (87) 


7.6 

8.0 

8.9 

10.0 

10.0 

10.8 

12.0 

12.5 

Stelwach  (88) 


.00187 

.00182 

.00180 

.00177 

.00171 

.00170 

.00169 

.00168 

.00187 

.00182 

.00180 

.00177 

.00171 

.00170 

.00169 

.00168 


.001955 

.001995 

.002117 

.002610 

.002559 

.001885 

.002095 

.002630 

.001767 

.001994 

.002575 

.001651 

•001872 

.001872 

.002494 

•001492 

.001802 

.002215 

.001495 

.001708 

.002160 

.001407 

.001594 

.001325 

•001452 

.002016 

.001243 

•001243 

.001293 

.002057 

.002057 
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M 


c 


1.897 

1.897 

1.897 

1.897 

1.897 

2.121 

2.121 

2.121 

2.121 

2.003 

2.003 

2.003 

2.003 

2.003 

2.186 

2.186 

2.186 

2.186 

2.186 


Rq  Rx  *  10”6  Cg  (exp)  Cd  (exp) 

Wilson  (75) 

3.1  .0033 

5.8  .0026 

H.l  .00238 

13.8  .00242 

16.6  .00230 

-8**33  .00350 

®*3  .00240 

11.0  .00220 

13.5  .00218 

2.83  .00310 

5*45  .00258 

18.6  .00234 

13.3  .00230 

18*0  .00220 

8.8  .00240 

18.5  .00219 

13.8  .00218 

15.4  .00212 

8.82  .00390 


Saltsnwn  end  Fisher  (60s) 


0.51 

24200 

0.81 

51900 

0.60 

23500 

0.90 

32100 

0.90 

21300 

1.26 

31300 

1.53 

27000 

1.36 

33600 

0.82 

13500 

0.64 

23200 

1.72 

60300 

0.81 

116100 

1.54 

106400 

31.1 

.00224 

49.4 

.00188 

30.0 

.00224 

31.5 

.00187 

22.3 

•002G4 

30.1 

.00200 

27.3 

.00195 

29.5 

.00199 

15.4 

.00226 

22.2 

.00218 

42.8 

.00149 

74.2 

.00171 

69.4 

.00142 
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APPENDIX  IV 


EXPERIMENTAL  DATA  OF  SKIN  FRICTION  COEFFICIENTS  AT  VARIOUS 
REYNOLDS  NUMBERS,  MACH  NUMBERS  AND  TEMPERATURE  RATIOS 
(For  Heat  Transfer) 


M 

e 

f  A 
w  e 

Ro 

R  X  10**6 

X 

Cf  (exp) 

Abbot 

(U 

3.9 

1 

5 

7.25 

1.8 

5 

7.25 

1.8 

7.5 

Hill 

(31) 

8.99 

7.68 

1245 

.000790 

9.04 

7.97 

1607 

.000891 

9.07 

8.28 

1908 

.000851 

9.10 

8.69 

2287 

.000800 

8.22 

7.17 

2081 

.000924 

8.25 

7.26 

2498 

.000910 

8.27 

7.34 

2885 

.000870 

8.29 

7.37 

3202 

.000820 

8.29 

7.41 

3451 

.000771 

■ill 

(32) 

8.25 

6.12 

2200 

3.3 

.000910 

8.27 

6.17 

2505 

3.7 

.000840 

8.28 

6.25 

2760 

4.2 

.000796 

8.29 

6.18 

2965 

4.7 

.000734 

9.04 

8.07 

1936 

3.1 

.000913 

9.07 

8.30 

2276 

3.P 

.000870 

9.10 

8.65 

2710 

4.7 

.000805 

10.03 

9.76 

1300 

l.C 

.000841 

10.04 

9.32 

1450 

2.!. 

.000761 

10.05 

8.91 

1680 

2 .( 

.000696 

10.06 

8.99 

1700 

2.(5 

.000673 

Hopkins,  Keener  and  Lottlo  (34) 


6.5 

2.843 

2262 

.00157 

6.5 

3.493 

2389 

.00152 

6.5 

2.692 

4&5 

.00122 

6.5 

3.601 

3300 

.00125 

(exp) 


.00239 

.00116 

.00104 
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M 

c 

t  /r 
w  e 

R-  R  ::  10~6 

6  x 

C£  (exp) 

Hopkins, 

Keener  and  Louie  (cont.) 

6.5 

3.503 

5900 

.00120 

6.5 

3.765 

3815 

.00123 

6.5 

4.330 

2185 

.00141 

6.5 

4.176 

3890 

.00135 

6.5 

4.300 

8326 

.00100 

6.5 

4.333 

6419 

.003 16 

Lee,  Yanta  and  Leonas  (42) 

4.72 

3.67 

38100 

.000770 

4.71 

3.65 

34450 

.000782 

4.71 

3.66 

31000 

.000792 

4.71 

3.66 

27400 

.000794 

4.70 

3.65 

23800 

.000812 

4.70 

3.64 

19750 

.000832 

4.695 

3.64 

15700 

.000856 

4.66 

3.58 

11340 

.000892 

4.59 

3.48 

6500 

.000936 

4.77 

2.78 

19180 

.000928 

4.79 

2.40 

17450 

.000976 

4.72 

3.69 

41750 

.000760 

Monaghan  and  Cooke  ( 5  o ) 


2.43 

2.94 

1893 

.98 

.00300 

2.43 

2.94 

2278 

1.40 

.00273 

2.43 

2.94 

3270 

2.10 

.00232 

2.43 

2.94 

3516 

2.51 

.00200 

2.43 

2.94 

4525 

3.20 

.00188 

2.43 

3.42 

2239 

1.42 

.00248 

2.43 

3.42 

3003 

2.12 

.00232 

2.43 

3.42 

3568 

2.71 

.00218 

2.43 

3.42 

4556 

3.70 

.00209 

2.43 

3.42 

5581 

4.39 

.00195 

Monaghan  and  Cooke (51a) 

2.82 

3.5 

1944 

1.343 

.00221 

2.82 

3.5 

2464 

1.731 

.00206 

2.82 

3.5 

3107 

2.426 

.00198 

2.82 

3.5 

3391 

2.897 

.00196 

Pappas 

(57) 

1.69 

1.7 

.807 

1.69 

1.7 

1‘.50 

1.69 

1.7 

2.20 

1.69 

1.7 

2.89 

1.69 

1.7 

3.55 

1.69 

1.65 

2.74 

Cd  (exp) 


.00388 

.00326 

.00313 

.00279 

.00281 

.00315 

.00280 

.00265 

.00245 

.00255 


.00290 

.0028!: 

.0025c 

.0023' 


.0036- 

.0032: 

.00316 

.0029: 

.00280 

.00310 
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H 

e 

T  /r 
w  e 

Rn  R  x  10“6 

0  x 

Cf  (exp) 

cd  (exp! 

Pappas  (cont.) 

1.69 

1.65 

3.80 

.00300 

1.69 

1.65 

6.00 

.002/6 

1.69 

1.63 

1.89 

.00336 

1.69 

1.63 

1.92 

.00337 

1.69 

1.63 

2.92 

.00310 

1.69 

1.63 

3.00 

.00316 

1.69 

1.63 

4.01 

.00283 

1.69 

1.63 

4.09 

.00271 

1.69 

1.63 

5.42 

.00279 

1.69 

1.61 

3.62 

.00273 

1.69 

1.61 

5.45 

.00255 

1.69 

1.61 

7.18 

.00256 

1.69 

1.61 

7.21 

.00240 

1.69 

1.61 

8.98 

.00235 

2.27 

2.18 

0.94 

.00342 

2.27 

2.18 

1.95 

.00307 

2.27 

2.18 

2.88 

.00286 

2.27 

2.18 

3.90 

.00270 

2.27 

2.18 

4.80 

.00270 

2.27 

2.19 

1.48 

.00340 

2.27 

2.19 

2.35 

.00290 

2.27 

2.19 

4.30 

.00268 

2.27 

2.19 

5.50 

.00254 

2.27 

2.12 

,  3.56 

.00273 

2.27 

2.12 

4.90 

.00250 

2.27 

2.12 

6.40 

.00250 

2.27 

2.12 

7.80 

.00230 

Sommer  and  Short  (65) 

2.81 

1.03 

3.00 

.00312 

3.82 

1.05 

4.07 

.00250 

5.63 

1.29 

4.71 

.00187 

6.90 

1.70 

4.06 

.00138 

6.90 

1.70 

6.09 

.00144 

7.0 

1.75 

6.06 

.00126 

7.0 

1.75 

9.92 

.00132 

3.78 

1.05 

4.94 

.00229 

3.67 

1.05 

3.78 

.00251 

Winkler  (76) 

5.26 

5.52 

3880  5.04 

.00135 

.00154 

5.29 

5.57 

4300  5.94 

.00131 

.00145 

4.98 

4.51 

1900  2.29 

.00134 

.00165 

5.19 

4.74 

1782  2.58 

.00161 

.00133 

5.20 

4.83 

2960  3.81 

.00135 

.00155 

5.24 

5.02 

3455  4.88 

.00115 

.00152 

5.24 

4.97 

3799  5.11 

.00106 

.00143 

5.17 

3.89 

1055  2.01 

.00147 

.00133 

5.16 

3.71 

1652  2.57 

.00132 

.00129 
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Rx  x  10 


6 


M 


c 


5.10 

5.11 
5.20 

,12 
,21 
14 
20 


Winkler  (cont.) 


3.58 

3.52 

3.77 

3.78 
5.15 
5.50 
5.38 


1735 

2488 

2482 

3256 

2099 

2936 

3173 


2.73 

2.27 

3.27 
3.57 
2.72 
3.36 
4.07 


Cf  (exp) 


.00134 

.00124 

.00120 

.00105 

.00147 

.00139 

.00143 


Cd  (exp) 


.00327 

.00153 

.00137 

.00182 

.00154 

.00175 

.00156 
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